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Scattered  sunlight  in  the  atmosphere  v/as  studied 
through  computer  modeling  and  comparison  with  aureole 
measurements.  An  analytical  technique  was  developed  to 
determine  the  aerosol  optical  depth  from  solar  radiometer 
measurements  in  the  ultraviolet.  A  photographic  technique 


was  used  to  measure  the  scattered  light  of  the  solar  aureole 
for  comparison  with  a  curved  Earth  modified  single  scat- 
tering calculation. 

Because  of  the  rapid  spectral  variation  of  the  ozone 
absorption  coefficients  below  350nm,  the  bandwidth  of  a 
filter  radiometer  (sunphotometer)  is  sufficient  to  inval- 
idate the  standard  Bouguer-Langley  analysis  used  to  deter- 
mine aerosol  optical  depth.  By  carefully  modeling  our 
instrument  with  O.lnm  wavelength  intervals,  we  were  able  to 
reproduce  the  observed  deviation  from  straight  line  Langley 
plots.  Using  this  model  we  found  an  analytic  fxinction  to 
characterize  the  deviation  from  linearity,  thus  allowing 
direct  computation  of  the  aerosol  optical  depth  even  in 
the  ultraviolet. 

A  modified  single  scattering  calculation  was  used  to 
predict  sky  brightness  near  the  sun  as  accurately  as  pos- 
sible without  going  to  the  formalism  of  a  full  multiple 
scattering  calculation.  The  calculation  uses  spherical 
Earth  geometry  but  does  not  include  refraction.  It  includes 
the  contribution  to  sky  brightness  from  sunlight  reflected 
from  the  ground  and  contains  an  empirical  correction  term 
for  Rayleigh  multiple  scattering.  Using  the  tables  of 
Coulson,  Dave,  and  Sekera  for  total  scattering  in  a  pure 
Rayleigh  atmosphere,  correction  terms  for  multiple  scat- 


vx 


tering  were  derived  for  the  Rayleigh  component  as  a  function 
of  wavelength  and  solar  zenith  angle.  Analytic  functions, 
were  used  to  represent  the  altitude  distribution  of  atmos- 
pheric constituents  based  on  the  measurements  of  Elterman. 

The  photographic  technique  for  aureole  study  is  de- 
scribed. The  sun  is  blocked  by  a  neutral  density  4.0 
occulting  disc,  allowing  the  sun  to  be  photographed  simul- 
taneously with  the  surrounding  sky  for  accurate  determina- 
tion of  the  sky  to  sun  intensity  ratio.  Film  sensitometry 
is  discussed,  including  the  use  of  a  new  analytic  function 
to  describe  the  characteristic  curves  of  film. 

A  series  of  graphs  of  calculated  and  measured  aureoles 
show  the  comparative  accuracy  and  range  of  validity  of  the 
modified  single  scattering  calculation.  The  calculation 
appears  to  accurately  describe  the  sky  brightness  except 
when  aerosol  multiple  scattering  becomes  dominant,  such  as 
at  middle  wavelengths  near  sunset.  At  short  wavelengths 
Rayleigh  scattering  is  dominant  and  the  Rayleigh  correction 
factor  aids  accuracy.  The  ultraviolet  (309nm)  aureole  is 
compared  with  a  Monte  Carlo  multiple  scattering  calculation, 
and  with  our  modified  single  scattering  calculation,  to 
show  that  the  latter  is  fairly  accurate  even  in  the  ultra- 
violet. 
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CHAPTER  1 


INTRODUCTION 


The  main  features  of  the  daytime  sky  brightness  were 
explained  long  ago  by  Rayleigh  as  scattering  of  sunlight 
by  molecules  (actually  by  molecular  density  fluctuations) 
in  the  atmosphere,  and  an  idealized  atmosphere  consisting 
only  of  air  molecules  is  known  as  a  Rayleigh  atmosphere. 
But  the  complexities  of  observed  sky  radiation  depend 
to  a  large  extent  on  the  atmospheric  aerosols,  the  partic- 
ulate matter  suspended  in  the  atmosphere  including  water 
droplets,  dust,  smoke,  and  other  air  pollutants  solid  or 
liquid.  In  addition  absorption  is  an  important  factor 
in  the  middle  ultraviolet  and  in  the  infrared,  but  in  the 
range  400nm  to  800nm  radiation  reaching  the  ground  under- 
goes extinction  primarily  by  scattering.  Thus  the  study 
of  light  scattering  by  atmospheric  aerosols  is  of  great 
interest  in  this  wavelength  range. 

Typical  aerosols  include  particles  with  radii  ranging 
from  lOnm  to  10,000nm,  and  although  the  range  of  sizes 
may  extend  well  beyond  these  limits ,  the  fact  that  the 
optical  wavelengths  occupy  the  range  400nm  to  700nni  makes 


these  particles  of  similar  size  important  in  light  scat- 
tering. An  important  parameter  used  widely  in  Mie  scat- 
tering calculations  is  the  size  parameter  x  =  2TTr/X 
where  r  is  the  particle  radius  and  X  is  the  wavelength  of 
the  scattering  radiation.  For  x  >  1.0  the  departure  from 
Rayleigh  scattering  is  generally  considerable. 

The  theory  of  scattering  of  electromagnetic  radiation 
from  spheres  of  arbitrary  size  and  index  of  refraction 
is  generally  credited  to  Mie  (1) ,  although  N.A.  Logan  (2) 
observes  that  Ludwig  Lorenz  published  a  similar  solution 
to  Mie's  eighteen  years  earlier  in  1890.  The  best  known 
book  on  Mie  scattering  is  one  by  Van  de  Hulst  (3) .  The 
Mie  solution  consists  of  a  series  whose  terms  contain 
Bessel  functions  of  half-integral  order  (spherical  Bessel 
functions)  with  complex  argument,  and  first  and  second 
derivatives  of  the  Legendre  polynomials.  The  number  of 
terms  required  for  reasonably  accurate  convergence  of  the 
series  is  on  the  order  of  the  size  parameter  x.  Only  since 
the  advent  of  the  high  speed  computer  has  it  been  possible 
to  perform  the  extensive  calculations  needed  to  describe 
the  wide  range  of  sizes  present  in  the  atmospheric  aerosol. 
The  Mie  scattering  calculations  for  this  dissertation  were 
done  using  a  subroutine  developed  by  Dave  (4)  which  uses 
a  downward  recurrence  procedure  that  accurately  gives 


Mie  scattering  intensities  even  for  very  large  particles 
of  x=200. 

An  additional  complexity  of  observed  sky  brightness 
is  multiple  scattering,  which  will  be  important  whenever 
the  intensity  of  scattered  radiation  is  a  non-negligible 
fraction  of  the  intensity  of  direct  sunlight  entering  a 
unit  voltome.  Such  a  condition  will  occur  when  the  aerosol 
optical  depth  is  large,  near  sunset  when  light  travels 
over  a  very  long  path  through  the  atmosphere,  and  for  the 
middle  ultraviolet  wavelengths  where  the  Rayleigh  optical 
depth  is  large-  Chandrasekhar  (5)  developed  the  theory 
of  multiple  scattering  for  a  Rayleigh  atmosphere  and  in 
1957  Sekera  (6)  used  this  theory  to  calculate  the  distri- 
bution of  Rayleigh  sky  radiation.  Application  of  multiple 
scattering  calculations  to  aerosols  is  more  difficult,  but 
Plass  and  Kattawar  (7)  have  used  a  Monte  Carlo  technique 
in  which  the  computer  does  a  light  scattering"experiment. " 
The  computer  injects  photons  into  the  atmosphere  and  fol- 
lows them  down,  deciding  when  scattering  will  occur  and 
what  kind  by  generating  a  random  number  and  applying  it  to 
a  weighted  probability  for  scattering,  non-scattering,  or 
absorption.  By  following  a  large  enough  niomber  of  photons 
for  good  statistics,  a  sky  brightness  distribution  results 
which  is  quite  accurate.  A  Monte  Carlo  program  based  on 


that  of  Plass  and  Kattawar  was  written  by  Furman  (8)  of 
our  group,  and  this  program  was  used  for  the  multiple  scat- 
tering results  in  this  work. 

The  composition  and  size  and  altitude  distribution 
of  the  atmospheric  particulate  matter  is  studied  by  a 
number  of  techniques.  A  good  survey  of  the  different 
methods  is  presented  in  a  book  by  Allen  (9)  with  most 
experimental  methods  falling  into  two  major  categories, 
direct  sampling  and  remote  sensing. 

Direct  sampling  involves  moving  a  collector  through 
the  volume  of  air  of  interest  and  analyzing  the  particles 
collected.  Ground  based  sampling  systems  can  study  only  the 
lowest  layer  of  the  atmosphere,  while  airborne  systems 
generally  cannot  sample  continuously  and  have  problems 
with  contamination  and  calibration.  It  is  interesting  to 
note  that  during  coordinated  studies  performed  at  the 
University  of  Florida  in  19  71  a  serious  discrepancy 
between  the  size  distribution  obtained  by  our  ground  based 
aureole  measurements  and  that  obtained  by  an  airborne 
particle  counter  was  resolved  with  recalibration  of  the 
particle  counter. 

Remote  sensing  determines  parameters  of  the  atmos- 
pheric aerosols  by  observing  the  effects  of  aerosols  on 
natural  or  manmade  radiation  fields.  This  dissertation 


concerns  measurements  of  natural  siinlight,  direct  and 
scattered  at  several  different  wavelengths,  from  which  we 
determine  the  aerosol  optical  depth,  number  density,  size 
distribution,  and  to  some  extent  composition  and  altitude 
distribution.  Experimental  work  concentrated  on  measure- 
ment of  the  attenuated  direct  solar  radiation  as  a  function 
of  wavelength  and  solar  zenith  angle  to  determine  aerosol 
optical  depths,  and  on  measurement  of  the  solar  aureole, 
the  region  of  intense  forward  scattered  light  in  the 
vicinity  of  the  sun,  to  determine  the  aerosol  size  distri- 
bution. 


CHAPTER  2 
AEROSOLS 

The  characteristics  of  atmospheric  aerosols  are 
important  because  they  affect  light  scattering  in  the 
atmosphere,  especially  the  scattered  light  of  the  solar 
aureole.  Unfortunately  aerosols  may  vary  in  a  number  of 
ways:  in  index  of  refraction,  chemical  composition,  size 
distribution,  and  shape  and  structure.  And  because  aerosol 
particles  are  constantly  being  produced  and  removed,  all 
these  parameters  may  vary  with  time.  In  order  to  draw 
conclusions  about  aerosol  particles  from  the  limited 
information  available  from  optical  study  of  the  solar 
aureole,  it  is  necessary  to  identify  the  important 
parameters  and  realistically  model  them  using  the  best 
information  available. 

2.1 Aerosol  Size  Distribution 

The  natural  aerosol  consists  of  all  solid  or  liquid 

colloidal  particles  present  in  the  atmosphere,  ranging  in 

size  from  the  molecular  complexes  of  the  small  and  inter- 


mediate  ions  to  the  giant  nucleii  20  microns  or  more  in 
radius.  (Here  I  shall  use  the  micron  as  the  unit  of  aerosol 
radius,  following  the  convention  of  aerosol  study.)  The 
radius  interval  of  greatest  interest  is  that  of  the  large 
nucleii,  from  0.1  micron  to  1.0  micron.  Particles  in  this 
size  range  are  most  effective  in  scattering  light. 

The  very  small  particles,  the  Aitken  nucleii,  are  not 
very  important  because  of  their  relatively  low  concen- 
tration and  because  they  scatter  light  ineffectively. 
The  theory  of  Smoluchowsky  describes  the  process  by  which 
small  particles  coagulate  to  form  large  particles.  Small 
particles  coagulate  more  quickly  than  large  because  of 
their  greater  nvimbers,  so  that,  for  instance,  0.01  micron 
particles  will  be  reduced  to  a  negligible  nvimber  in  only 
half  a  day. 

For  particles  larger  than  the  Aitken  nucleii,  the  num- 
ber of  particles  per  unit  volume  in  the  size  range  r  to 
r+dr  generally  decreases  with  radius  as  r~^  where  the 
exponent  of  the  decrease  v  is  found  to  be  1.3  as  a  world- 
wide average.  Such  a  distribution  ia  known  as  a  Junge  (10) 
distribution. 

Sedimentation,  direct  fallout,  and  perhaps  washout 
by  rain  decrease  the  number  of  giant  nucleii  to  the  point 
that  they  are  not  important  in  atmospheric  light  scattering. 
Thus,  to  a  good  approximation,  scattering  calculations  may 


be  performed  using  a  simple  Jiinge  distribution  with  an 
imposed  lower  limit  of  0.04  micron  or  so  and  an  imposed 
upper  limit  of  10.0  microns. 

Many  aerosol  particles  are  produced  through  conden- 
sation. According  to  Bullrich  (11)  ,  combustion  processes 
are  responsible  for  the  formation  of  the  very  small 
particles,  with  the  gaseous  combustion  products  condensing 
upon  cooling  to  form  smoke.  Subsequent  coagulation  then 
produces  the  larger  atmospheric  particles,  which  are  mainly 
man  made.  Particles  of  condensation  origin  may  be  either 
solid  or  liquid.  Hygroscopic  materials  such  as  ammonium 
chloride  may  undergo  a  phase  transition  from  solid  to 
liquid,  depending  on  the  humidity  of  the  medium. 

Other  aerosol  particles  are  produced  by  dispersion, 
the  process  of  atomization  of  solids  and  liquids  and  their 
injection  into  the  atmosphere  by  the  action  of  winds  or 
vibration.  Dispersion  aerosols  are  generally  solid  and  of 
varied  shape. 

Specific  mathematical  models  for  the  aerosol  size 
distribution  will  be  discussed  in  Chapter  5. 


2.2 Physical  Properties  of  Aerosols 

Inasmuch  as  Mie  scattering  calculations  assume  that 
all  particles  are  spherical,  the  effect  of  irregular 
aerosol  particles  must  be  considered.  As  mentioned  above, 
the  larger  aerosol  particles  produced  by  dispersion  are 
more  likely  than  the  small  condensation  particles  to  be  of 
irregular  shape.  In  addition,  ice  crystals  are  known  to 
form  as  hexagonal  rods  and  platelets,  and  salt  crystals 
formed  by  evaporation  of  sea  spray  are  cubic.  The  conclu- 
sions drawn  by  Holland  and  Gagne  (12) ,  that  scattering 
from  irregular  particles  deviates  from  that  of  spherical 
particles  only  for  large  particles  and  only  for  large 
scattering  angles,  indicate  that  particle  irregularity 
will  not  be  an  important  factor  for  this  dissertation, 
which  concentrates  on  light  scattered  at  small  angles. 
To  quote  Holland  and  Gagne,  "...the  unpolarized  mass  scat- 
tering coefficient  for  a  polydisperse  system  of  irregular 
but  randomly  oriented  particles  shows  a  remarkable  similar- 
ity to  the  corresponding  coefficient  for  spherical  parti- 
cles, but  ...  predicting  backscattering  for  such  particles 
on  the  basis  of  the  spherical  model  can  lead  to  serious 
errors"  (p.  1120) . 

The  index  of  refraction  is  a  further  parameter  that 
affects  light  scattering.  For  years  an  index  of  refraction 
of  1.33,  that  of  water,  was  assumed  for  lack  of  any  better 
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estimate;  but  more  recent  data  indicate  that  an  index  of 
refraction  of  1.50  to  1.55  is  more  accurate.  Volz  (13) 
has  estimated  the  index  of  refraction  from  distillation 
residue  of  rainwater  and  obtained  a  value  of  1.53,  but 
this  method  is  applicable  only  to  the  water  soluble 
component,  which  typically  constitutes  only  about  15% 
of  the  total  aerosol.  Estimates  have  also  been  made  based 
on  cuialysis  of  chemical  composition  of  collected  aerosol. 
Recently  Bhardwaja  (14)  has  used  integrating  nephalometer 
measurements  of  total  scatter  and  hemispheric  backscatter 
for  an  in  situ  measurement  of  the  refractive  index,  and 
obtained  an  index  of  refraction  of  1.55  ±  0.03  with  an 
imaginary  component  of  0.03. 

Another  method  for  determining  the  index  of  refraction 
is  use  of  bistatic  laser  scattering  to  measure  the  polar- 
ization of  light  over  a  large  range  of  scattering  angles. 
This  technique  was  used  by  G.  Ward  (15)  at  the  University 
of  Florida  in  19  72,  who  obtained  an  index  of  refraction 
of  1.50  +  O.OOSi. 

Because  of  the  insensitivity  of  forward  scattering  to 
index  of  refraction,  no  attempt  should  be  made  to  infer 
it  from  study  of  the  solar  aureole.  On  the  other  hand, 
this  insensitivity  means  that  assumption  of  a  reasonable 
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value  for  the  index  of  refraction,  1.50+O.Oli,  will  not 
lead  to  serious  error  in  aureole  calculations. 

Parameters  that  affect  the  index  of  refraction  of 
aerosols  are  the  chemical  composition  and  the  humidity.  As 
noted  earlier,  condensation  aerosols  are  frequently  hygro- 
scopic and  undergo  a  phase  transition  from  solid  to  liquid 
with  increasing  humidity.  According  to  data  presented  by 
Bullrich  (11)  such  a  particle  with  an  index  of  refraction 
of  1.5  will  show  little  change  until  a  humidity  of  80% 
is  reached,  but  between  80%  and  100%  humidity  the  index 
of  refraction  rapidly  decreases  to  that  of  water,  1.33. 
Simultaneously  the  particle  grows  in  size,  but  this  is  a 
smaller  effect.  At  100%  humidity  a  typical  hygroscopic 
particle  is  only  1.3  times  larger  than  in  its  solid  phase. 

2.3 Aerosol  Altitude  Distribution 

The  aerosol  number  density  is  found  to  decrease  very 
rapidly  with  altitude,  an  order  of  magnitude  in  the  first 
three  kilometers,  but  is  found  also  to  have  a  high  altitude 
layer  at  approximately  20  kilometers.  The  bulk  of  all  aero- 
sols are  found  in  the  first  few  kilometers  and  they  are 
generally  as  described  in  the  preceding  sections.  The 
microstructure  of  the  tropospheric  aerosol  distribution 
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has  been  studied  by  Rosen  (16)  who  found  that  the  altitude 
distribution  can  be  far  more  complicated  than  a  simple 
exponential  falloff ,  especially  when  temperature  inversions 
are  present.  Layering  and  variation  in  size  distribution 
as  well  as  number  density  were  found.  In  some  cases  the 
mixing  layer  model,  in  which  particles  are  assumed  vuiiform- 
ly  distributed  up  to  some  mixing  height,  may  be  an  accurate 
model.  To  a  large  extent,  however,  the  fact  -that  light 
from  the  sun  must  traverse  the  entire  atmosphere  will 
average  out  the  effects  of  these  small  scale  variations. 

The  stratospheric  dust  layer  has  been  studied  by 
Junge  (17)  using  balloon  and  aircraft  borne  detectors. 
He  found  that  particles  less  than  0.1  micron  decreased 
rapidly  in  number  with  height  and  were  of  tropospheric 
origin.  The  peak  of  the  particle  size  distribution  was 
between  0.1  and  1.0  microns  with  a  number  density  of  about 
1.0  per  cubic  centimeter.  Chemical  analysis  showed  that 
they  were  mostly  sulphates  of  ammonia  and  were  probably 
formed  in  the  stratosphere  by  oxidation  of  gaseous  sulphur 
compounds.  Particles  larger  than  1.0  microns  were  very 
rare  and  were  distributed  according  to  a  different  power 
law,  suggesting  that  they  were  of  extraterrestrial  origin. 
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Volz  (18)  has  used  twilight  measurements  to  make  long 
terra  studies  of  the  stratospheric  aerosol  layer.  As  the 
sun  sets  the  shadow  of  the  edge  of  the  Earth  sweeps 
upward  through  the  atmosphere , , and  the  presence  of  a  dust 
layer  will  be  revealed  as  a  sudden  change  in  sky  brightness. 
His  studies  show  the  importance  of  volcanic  eruptions  in 
injecting  large  quantities  of  dust  into  the  stratosphere, 
dust  that  may  persist  for  more  than  a  year. 

Several  workers  have  used  active  source  remote  sensing 
systems,  high  intensity  laser  beams  or  searchlight  beams, 
to  study  the  aerosol  altitude  distribution.  Elterman  (19) 
directs  a  searchlight  beam  vertically  into  the  atmosphere 
and  observes  the  scattered  light  emanating  from  the  common 
volume  formed  by  the  intersection  of  the  searchlight  beam 
and  the  observation  cone.  By  scanning  along  the  beam  he 
obtains  information  on  aerosol  density  at  progressively 
greater  altitudes,  as  high  as  50  kilometers.  The  altitude 
distributions  of  air,  ozone,  and  aerosol  used  for  compu- 
tations in  this  dissertation  are  based  on  the  standard 
atmosphere  tabulated  from  the  measurements  of  Elterman. 


CHAPTER  3 
SOLAR  RADIATION  IN  THE  ATMOSPHERE 

3.1     Units 

In  view  of  the  confusion  that  has  traditionally 
existed  in  the  fields  of  photometry  and  radiometry,  a 
review  of  radiometric  units  currently  being  used  might  be 
useful  at  this  point. 

The  confusion,  inconsistency,  and  duplication  in  the 
definitions  of  fundamental  quantities  arise  in  part  from 
the  fact  that  the  fields  of  radiometry  and  photometry 
were  for  years  independent  and  not  clearly  connected. 
Radiometry  deals  with  measurement  of  electromagnetic 
energy  regardless  of  wavelength,  v/hile  photometry  is 
restricted  to  visible  light  wavelengths.  Conversion  of 
radiometric  values  to  photometric  must  take  into  account 
the  response  of  the  human  eye.  This  response,  called  the 
luminous  efficiency  curve,  is  set  to  1.0  at  555  nanometers 
and  decreases  to  0.001  at  410nm  and  720nm.  The  conversion 
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factor  between  photometric  and  radiometric  units  is  that 
one  watt  of  monochromatic  555nm  light  is  equivalent  to 
685  lumens.  Measurements  at  other  wavelengths  must  be 
weighted  by  the  luminous  efficiency  curve . 

The  radiant  quantities  which  are  employed  in  this 
dissertation  follow  a  summary  of  basic  radiometric  quanti- 
ties defined  by  McCluney  (20)  in  1968.  These  quantities 
are  summarized  in  Table  3.1. 

The  radiant  flux  may  be  defined  as  the  radiant  energy 
per  unit  time  and  is  measured  in  watts.  The  radiant  energy 
per  unit  time  per  unit  wavelength  interval  is  called  the 
spectral  radiant  flux  and  is  measured  in  watts  per  nano- 
meter. The  term  "spectral"  will  always  denote  per-unit- 
wave length  measurements.  The  photometric  quantity  corres- 
ponding to  the  radiant  flux  is  the  luminous  flux,  which 
is  measured  in  lumens.  If  P  is  the  spectral  radiant  flux 
and  K  is  the  luminous  efficiency  curve,  the  luminous  flux 
F  in  lumens  is  given  by 

F  =  685  /"  K(A)  P(A)  dA  (3.1) 

Intensity  is  a  term  that  is  often  misused  and  over- 
worked. The  radiant  intensity  is  defined  as  the  radiant 
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TABLE  3.1 


Radiometric  Units 


■1 


Name  (spectral-) 

Radiant  flux 
Radiant  intensity 
Radiance 
Irradiance 


Units 


watts/nm 

watts/ (nm  sr) 

2 

watts/ (nm  sr  cm  ) 

2 

watts/ (nm  sr  cm  ) 


Photometric  Units 


Luminous  flux 

lumens 

Luminous  intensity 

candle=lumen/sr 

Luminance 

2 
candle/cm 

Illuminance 

2 

lumen/cm 
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flux  per  unit  solid  angle  from  a  point  source.  A  similar 
term  applicable  to  extended  sources  is  the  radiance,  which 
is  radiant  flux  per  unit  solid  angle  per  unit  area  of 
emitting  surface.  Corresponding  photometric  terms  are 
Iximinous  intensity  and  luminance,  measured  in  terms  of  the 
candle,  which  is  one  lumen  per  steradian,  and  candles  per 
unit  area  respectively.  These  are  calculated  in  a  similar 
fashion  to  equation  3.1. 

Finally,  the  irradiance  is  defined  as  the  radiant  flux 
incident  upon  a  surface  per  unit  area.  The  magnitude  of 
irradiance  from  a  point  source  follows  an  inverse  square 
decrease  with  distance;  the  irradiance  from  an  extended 
source  may  in  some  cases  decrease  more  nearly  linearly. 
The  cosine  law  states  that  if  either  an  emitting  or  an 
irradiated  surface  is  other  than  perpendicular  to  the 
direction  of  the  light,  the  irradiance  will  vary  as  the 
cosine  of  the  angle.  Irradiance  is  measured  in  watts  per 
square  centimeter  and  the  corresponding  photometric  unit, 
the  illuminance,  in  limens  per  steradian  per  unit  area. 

By  these  definitions  a  sunphotometer,  which  measures 
the  intensity  of  sunlight  over  a  narrow  wavelength  band 
after  transmission  through  the  atmosphere,  should  more 
properly  be  called  a  solar  radiometer,  especially  when 
the  instrument  operates  at  non-visible  wavelengths. 
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3.2 The  Extraterrestrial  Solar  Flux 

Since  this  work  concerns  remote  sensing  of  the  atmos- 
phere using  the  sun  as  the  light  source,  it  is  appropriate 
at  this  point  to  discuss  some  of  the  features  of  this  light 
source. 

Viewed  at  a  distance  of  one  astronomical  unit,  the 
average  distance  from  the  Earth  to  the  sun,   the  sun  sub- 
tends thirty-two  minutes  of  arc,  little  enough  for  the 
finite  extent  of  the  source  to  be  neglected  for  many  radio- 
metric calculations.  Light  intensity  is  fairly  constant 
across  the  face  of  the  sun  in  this  wavelength  region, 
varying  about  30%  due  to  limb  darkening  and  with  some 
variation  due  to  sunspots,  small  scale  granularity,  and 
possible  solar  flares. 

Because  the  Earth's  orbit  is  elliptical,  there  will 
be  a  variation  in  intensity  amounting  to  about  3.4%  at 
most  because  of  the  varying  sun-Earth  distance,  with  the 
Earth  being  closest  to  the  sun  in  December. 

With  the  exception  of  the  very  short  wavelengths,  the 
energy  output  of  the  sun  is  quite  stable.  The  total 
irradiance  of  the  sun,  the  solar  constant,  was  measured 
by  Moon  (21)  in  1940  to  be  1323  watts  per  square  meter  and 
by  others  since  to  be  1374  in  1950,  1395  in  1954,  and  1361 
in  1968.  These  uncertainties  are  due  mostly  to  difficulties 
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in  making  measurements  through  the  highly  variable  dust, 
haze,  cuid  cloud  cover  of  the  Earth's  atmosphere.  In  1969 
Thekaekara  et  al.  (22)  reported  the  results  of  a  series 
of  solar  constant  measurements  made  from  an  aircraft  at  an 
altitude  of  11.6  kilometers,  above  80%  of  the  permanent 
gases  of  the  atmosphere  and  above  99.9%  of  the  highly 
variable  water  vapor,  dust,  and  smoke.  They  report  an 
average  solar  constant  of  1351  watts  per  square  meter. 

The  spectral  distribution  of  sunlight  may  be  approx- 
imated in  the  visible  by  a  smooth  7200 "K  continuum  with 
superimposed  Fraunhofer  absorption  lines,  and  as  a  6000 °K 
greybody  for  wavelengths  beyond  1250nm.  On  a  microscopic 
scale  siinlight  below  600nm  is  extremely  complicated,  so 
that  determining  the  continuum  is  difficult.  This  is 
especially  a  problem  between  300nm  and  450nm  and  is  one 
source  of  the  greater  uncertainty  in  the  solar  flux  in 
this  region. 

Probably  the  most  accurate  determinations  of  the 
extraterrestrial  spectral  solar  irradiance  were  made  on 
the  same  1968  aircraft  flights  that  the  solar  constant 
measurements  mentioned  above  were  made.  Although  Thekaekara 
(22)  and  Arvesen  et  al.  (23)  made  irradiance  measurements 
from  the  same  airplane  at  the  same  time,  Thekaekara  using 
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TABLE  3.2 


Solar  spectral  irradiance  at  selected  wavelengths 


A(nm) 

F(X)  w/m^  nm 

(smoothed) 

250 

0.068 

0.068 

300 

0.545 

0.560 

350 

1.015 

1.040 

400 

1.692 

1.540 

450 

2.023 

2.060 

500 

1.927 

1.960 

550 

1.849 

1.880 

600 

1.816 

1.790 

650 

1.577 

1.610 

700 

1.495 

1.500 

800 

1.154 

1.160 

900 

0.889 

0.893 

1000 

0.741 

0.741 
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a  Perkin-Elmer  monochroinator  and  Arvesen  et  al.  using  a 
Gary  14  monochroinator,  there  are  discrepancies  between  the 
two  solar  spectra.  Integrating  over  his  spectral  measure- 
ments Thekeakara  finds  a  solar  constant  of  1351  watts/m 

2 
while  Arvesen  et  al.  find  1390  watts/m  .  We  have  chosen  to 

use  the  spectrum  of  Arvesen  et  al.  because  that  of  Thek- 

aekara  fails  to  give  the  correct  wavelengths  for  several 

well  known  Fraunhofer  dips,  while  that  of  Arvesen  et  al. 

correctly  gives  the  Fraunhofer  structure.  No  explanation 

for  the  discrepancy  is  given,  but  it  might  be  caused  by 

scanning  with  a  large  slit  width. 

We  calculate  the  apparant  position  of  the  sun  in  the 

sky  using  ephemeris  data,  the  exact  latitude  and  longitude 

of  the  observation  site,  and  the  local  time.  In  Figure  3.1 

we  show  the  time  at  which  the  sun  rises,  crosses  the 

meridian,  and  sets  throughout  the  year  in  Gainesville, 

neglecting  refraction.  Also  shown  are  the  azimuth  angle 

of  the  sun  at  sunrise  and  the  noon  solar  zenith  angle. 
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3.3    Geometry  of  the  Atmosphere 

An  assumption  common  to  atmospheric  optics  calcula- 
tions is  that  of  horizontal  homogeneity,  so  that  optical 
properties  of  the  atmosphere  depend  on  only  one  coordinate, 
the  height  above  ground  level.  This  assumption  is  justified 
for  the  molecular  component  of  the  atmosphere,  but  less  so 
for  the  aerosol  component.  As  point  sampling  by  particle 
counters  shows,  there  can  be  large  variations  in  aerosol 
content  from  place  to  place  depending  on  local  sources. 
The  fact  that  optical  measurements  are  made  over  long  path- 
lengths  somewhat  alleviates  problems  of  horizontal  homogen- 
iety  through  averaging,  but  the  experimentalist  must  be 
aware  of  the  problem. 

The  atmospheric  air  mass  is  a  measure  of  the  amount 
of  air  traversed  by  sunlight  before  reaching  the  ground. 
If  light  propogating  vertically  through  the  -  atmosphere  to 
sea  level  is  defined  to  have  penetrated  one  air  mass,  then 
to  a  first  approximation  light  propogating  at  zenith  angle 
9  traverses  an  air  mass  of  m(8)  =  sec{e).  This  is  a  result 
of  the  flat  Earth  approximation,  which,  because  the  atmos- 
phere is  a  thin  layer  on  a  very  large  sphere,  gives  very 
accurate  air  mass  values  for  zenith  angles  from  zero  to 
about  seventy-five  degrees.  As  the  zenith  angle  approach- 
es ninty  degrees  the  flat  Earth  sec (6)  approximation  goes 
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/////  /////////// 


(a)  Flat  Earth  approximation. 


(b)  Spherical  geometry. 


Figure  3.2 
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to  infinity,  while  actually  due  to  curvature  the  air  mass 
only  increases  to  thirty  or  so.  Spherical  Earth  corrections 
for  the  air  mass  for  an  exponential  atmosphere  were  calcu- 
lated by  two  men,  Bemporad,  whose  work  was  tabulated  by 
Schoenberg  (24)  ,  and  Chapman  (25)  .  Rozenberg  (26)  has 
proposed  an  empirical  expression  for  the  air  mass. 

m(6)  =  l/(cos(0)  +  0.025  e"^^°°^^®^)         (3.2) 

Bemporad 's  results  along  with  Rosenberg's  empirical  formula 
are  compared  with  sec (6)  in  Table  3.3. 

The  air  mass  is  defined  as  the  ratio 

m(e)  =  />^^^  ^^  (3.3) 

/"'N(h)  dh 

where  N(h)  represents  the  height  distribution  of  absorbing 
or  scattering  material,  dh  is  a  differential  element  of 
vertical  path  length,  and  ds  is  a  differential  element  of 
actual  path  length. 

Figure  3.2b  shows  the  geometry  of  an  incoming  ray 
traversing  a  spherical  atmosphere  for  which 
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TABLE  3.3 


Air  mass  as  a  function  of  solar  zenith  angle. 


9   sec 8   Bemporad  Rozenberg  seq  6  seq  6  seq  6  seq  C 


0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

30 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

60 

2.00 

2.00 

2.00 

2.00 

2.00 

1.98 

2.00 

70 

2.92 

2.90 

2.92 

2.91 

2.91 

2.84 

2.90 

80 

5.76 

5.60 

5.65 

5.61 

5.66 

5.14 

5.60 

85 

11.5 

10.4 

10.4 

10.4 

10.7 

8.04 

10.3 

88 

28,7 

19.8 

19.4 

18.7 

20.8 

10.5 

18.2 

90 

oo 

35-40 

40 

24.3 

30.2 

11.2 

23.6 
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ds  =  dh  sec (6) 


(3.4) 


Snell's  law  written  in  spherical  coordinates  appears  as 


sin((l))  = 


r  n(0)  sin(e) 
(r+h)  n(h) 


(3.5) 


where  n(h)    is  the   index  of  refraction  as   a   function  of 
height,    r  is  the  radius   of  the  Earth,   h   is   the   altitude, 
(j)  is  the   zenith  angle  of  the  ray  at  height  h,    and   6   is  the 
zenith  angle  of  the   ray  at  the   ground.   According   to 
Kondratyev   (27)    a  first  order  expression   for  the   index  of 
refraction  is 


n(0)      _ 
HThT    - 


1  +  2a 


1  +  2a{p(0)/p(h)  } 


(3.6) 


where  a=2.9xio'*.  We  may  now  write 


ds  = 


dh 


(  1  -  (-^)^x(!li°)^)^.si„Me) 


(3.7) 
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We  may  now  use  this  expression  in  eqn.  3.3  to  give  a 
general  formula  for  the  air  mass. 


B(h)  dh 
m(8)  =  /  _ _  (3^g) 


^'-   i-^)'<^^)'^si.Ue)    }V2 


r+h  ^   »•  n(h) 
By  the  usual  definition 

T  =  r  B(h)  dh  (3.9) 


0 


where  the  altitude  distribution,  the  volume  scattering 
coefficients  3(h)  give  the  attenuation  per  kilometer  for 
a  given  altitude.  Analytic  expressions  for  the  3(h)  will 
be  given  in  the  next  chapter.  Using  those  distributions  we 
calculated  the  total  air  mass  and  also  the  individual  air 
mass  contributions  for  the  Rayleigh,  Mie,  and  ozone  compo- 
nents. Where  Bemporad  and  Chapman  assumed  an  exponential 
distribution,  our  analytic  distributions  are  based  on  the 
measurements  of  Elterman  and  accurately  model  the  atmos- 
phere . 

Using  our  calculated  air  masses,  we  represented  each 
of  the  constituent  air  mass  functions  by  an  analytic 
modification  of  the  secant  law  introduced  by  Green  (28) 
and  called  a  seq  function,  defined  as 
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seq.(e)  =  T  ._         (3.10) 

^      [  1  -  sinMe)/q.)i/2 


where    q^  =  1.0018     Rayleigh 
q2  =  1.0074     Mie 
q,  =  1.0003     ozone. 

Attenuation  through  the  atmosphere  is  thus  accurately 
given  for  the  three  species  by 

I  =  Iq  expi'l    T^  seq^(e)  )  (3.11) 

i 

The  individual  and  total  air  masses  for  the  atmospheric 
species  are  compared  to  secant  and  Bemporad  air  mass 
results  in  Table  3.3. 


CHAPTER  4 


SOLAR  RADIOMETRY 


The  techniques  usually  used  in  the  analysis  of 
sunphotometry  and  radioinetry  data  are  not  applicable  below 
a  wavelength  of  350nm  because  of  the  spectral  structure 
in  the  ozone  absorption  coefficient.  We  have  developed  (29) 
experimental  and  analytical  techniques  for  extending  the 
range  of  solar  radioinetry  to  the  middle  ultraviolet. 

4.1    The  Attenuating  Components 

As  sunlight  traverses  the  atmosphere  three  species 
remove  energy  from  the  direct  beam:  air  molecules  Rayleigh 
scatter  light,  aerosols  remove  light  by  Mie  scattering, 
and  gases  such  as  ozone  absorb  light.  In  the  infrared  other 
gases  such  as  water  vapor  and  carbon  dioxide  replace 
ozone  as  the  absorbing  component.  The  total  optical  depth 
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is  the  sum  of  the  individual  optical  depths . 


^R  -^  ^M  -^  -"oz  (4-1^ 


Because  the  optical  depth  due  to  aerosols  is  highly  variable 
and  has  a  large  effect  on  visibility,  it  is  the  term  of 
greatest  interest  in  solar  radiometer  studies.  The  Rayleigh 
and  ozone  terms  are  much  less  variable  and  are  subtracted 
from  the  total  optical  depth  to  obtain  the  aerosol  optical 
depth . 

4.1.1    Rayleigh  Scattering 

Rayleigh  scattering  is  the  scattering  of  electromag- 
netic radiation  by  particles  very  much  smaller  than  the 
wavelength  of  the  incident  radiation,  i.e.  the  scattering 
of  light  by  air  molecules .  The  intensity  of  radiation 
scattered  from  a  volume  of  air  molecules  depends  on  the 
intensity,  polarization,  and  wavelength  of  the  incident 
radiation,  on  the  molecular  number  density  and  index  of 
refraction  of  air  for  the  particular  wavelength,  and  on 
the  angle  between  the  incident  and  scattered  radiation. 
The  equation  for  Rayleigh  scattering  is  written 
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R(4',X)  =   ^  ^    (m  -  1)^  (1  +  cos^dp))         (4.2) 
No  ^' 

where  R  is  the  Rayleigh  scattering  coefficient  for  angle  \l> 
and  wavelength  \.    Nq  is  the  molecular  number  density  and 
m  is  the  index  of  refraction.  If  for  a  given  wavelength 
the  intensity  of  an  unpolarized  incident  beam  is  I_ ,  the 
scattered  intensity  is 


Ki/;)  =  R(»p,A)  Iq  (4.3) 

Equation  4.2  displays  the  l/X**  dependance  that  makes 
Rayleigh  scattering  very  important  in  the  blue  and  ultra- 
violet. The  strong  scattering  of  blue  light  out  of  direct 
sunlight  produces  the  diffuse  "blue  sky"  background. 

Rayleigh  scattering  does  not  precisely  obey  a  l/A**  law 
because  it  also  depends  on  the  index  of  refraction  of  air, 
which  is  also  wavelength  dependent.  According  to  Edlen  (30) 
the  index  of  refraction  of  air  at  standard  temperature  and 
pressure  is  given  by 


m  -  1  =  0.000064328  +  O-O^^^^^J-  +  O'OOO^^^^O         (4.4) 

146  -  l/X^    41  -  l/A^ 
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where  A  must  be  in  microns. 

The  Rayleigh  attenuation  coefficient  is  the  product 
of  the  molecular  optical  cross  section  and  the  number  of 
molecules. 

3j^(h)  =  Oj^(A)   Nj^(h)  (4.5) 

The  cross  section  is  given  by 

8  tt'  (m^  -  1)2    6  +  36 

a  (X)  = (4.6) 

^      3  X"  Nj         6-76 

where         «•!•,-  ^  •    ,   2 » 

o-  =  Rayleigh  cross  section  (cm  ) 

X  =  wavelength  (cm) 

m  =  index  of  refraction  of  air 

Nq  =  molecular  number  density  for  a 

-3 
standard  atmosphere  (cm  ) 

19 
Nq  is  usually  given  as  2.547x10   particles  per  cubic 

centimeter. 

The  altitude  distribution  of  air  is  approximately 

exponential,  but  it  is  more  accurately  represented  by  a 

form  of  the  generalized  distribution  function  (31) 
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(1  312)2   h/6.42 
^       "   (0.312  +  eV6.42)^ 


where  h  is  the  altitude  in  kilometers.  A  comparison  of  the 
altitude  distribution  calculated  using  eqns.  4.5  and  4.7 
with  the  altitude  distribution  given  by  Eltenran  is  shown 
in  Figure  4.1. 

4.1.2    Mie  Scattering 

Mie  scattering  is  the  scattering  of  electromagnetic 
radiation  by  spherical  particles  of  arbitrary  size; 
Rayleigh  scattering  is  simply  the  very-small-particle  limit 
of  Mie  scattering.  Recently  progress  has  been  made  on 
calculating  the  scattering  from  arbitrarily  shaped  parti- 
cles, but  since,  as  explained  earlier,  this  work  is 
concerned  with  small  angle  scattering,  the  classic  Mie 
solution  for  spherical  particle  scattering  is  adequate. 

According  to  the  Mie  solution  the  scattering  from  a 
particle  is  a  function  of  the  size  of  the  particle  relative 
to  the  wavelength  of  the  incident  light,  so  the  scattering 
depends  on  the  size  parameter  x  =  27rr/A.  A  good  treatment 
of  Mie  scattering  is  given  by  Van  de  Hulst  (3)  whose 
notation  we  follow. 

The  problem  is  formulated  as  a  plane  electromag- 
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Altitude  distributions  of  the  SOOnm  Rayleigh, 
aerosol,  and  ozone  attenuation  coefficients 
according  to  Elterman  (1968)  compared  with 
analytic  functions. 


Figure  4 . 1 
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netic  wave  incident  on  a  dielectric  sphere  of  index  of 
refraction  m.  Beginning  with  Maxwell's  equations,  the  form 
of  the  solution  is  postulated  in  terms  of  spherical  Bessel 
fxinctions  for  the  interior  and  exterior  of  the  sphere  with 
boundary  conditions  applied  at  the  surface.  The  solution 
is  in  terms  of  two  scalar  constants  composed  of  Ricatti- 
Bessel  functions. 


a  = 

n 


(4.8) 


j^  _  m  ipj^Cmx)  'l'j^(x)  -  ij'j^(mx)  ^^ix) 


n 


m  ^^(mx)    Cj^(x)  -  rii^(mx)    C^(x) 


(4.9) 


where  primes  indicate  differentiation  with  respect  to  r 
the  particle  radius.  The  equation  for  the  electromagnetic 
wave  scattered  from  the  particle  is 


El 
E2 


ik(z-r) 


ikr 


0 


0 
S. 


El, 


E2, 


(4.10) 


where  El  and  E2  are  the  components  of  the  electric  field 
parallel  and  perpendicular  respectively  to  the  scattering 
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plane.  The  S,  and  S2  terms  are  expressed  in  terms  of  the 

parameters  a„  and  b  . 
n      n 

"   -)„^i        P^(cose)      d   , 

S  =  I     .Jn+J^   (  a„  -^ -^  K  ~   P;(cose))   (4.11) 

^     n=l  n(n+l)     ^    sin8      "  d8  ^ 


*   '>r^4.^  Pi(cose)      d   , 

So  =  I   -^^^^^   (  l^n  -^^ +a„—  Pi(cos0)    (4.12) 

^     n=l  n(n+l)     ^  sine       "  de  ^ 


The  differential  scattering  cross  section  is 


"e<-'-.«  =-^  (  isj'   |s,|M       ,,.,3, 


When  multiplied  by  the  irradiance  incident  on  a  particle, 
Oq   gives  the  intensity  of  light  scattered  in  the  direction 
6.  The  units  of  Oq   as  defined  are  area  per  unit  angle. 
The  total  cross  section,  which  is  a  measure  of  the  total 
energy  which  will  be  removed  from  the  incident  beam  by 
one  particle, is  found  by  integrating  over  all  angles. 


a(x,m)  =  /  aQ(x,m,0)  dw  (4.14) 

4Tr   ^ 
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In  order  to  apply  Mie  theory  to  atmospheric  scattering 
work  the  scattering  for  individual  particles  must  be 
integrated  over  all  particles  found  in  the  atmosphere.  It 
would  be  more  accurate  to  allow  the  size  distribution  to 
vary  with  altitude,  but  we  simplify  the  problem  by  assuming 
that  the  size-altitude  distribution  is  separable.  The  Mie 
extinction  coefficient  is  found  by  integrating  over  the 
size  distribution. 

e(A,m,h)  =  /"  a(x,m)  Njj(h)  n(r)  dr      (4.15) 

For  our  size  distribution  n(r),  the  fraction  of 
particles  in  the  radius  range  r  to  r+dr,  we  use  a  regular- 
ized power  law  distribution  which  corresponds  to  the  Junge 
distribution  for  large  r. 

r^-1 
n(r)  =  —   (4.16) 

^    {  1  +  (r/a)''  }^ 

The  parameter  v  is  the  power  of  the  number  decrease  and 
a  is  approximately  the  radius  of  maximum  frequency. 

The  altitude  distribution  we  represent  as  the  sum  of 
two  generalized  distribution  functions,  the  second  term 
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being  included  to  represent  the  high  altitude  aerosol 
component. 


1.135  e^/1-18 

^M^h)  =  N  (0)   [  uTT-To 

^  M     I  ^  Q  Qgg^  ^  3h/1.18j2 


[0)  f 


8.347  e^/-^*-'     ) 
+  (4.17) 

(80.0  +  e^/^-^)2 


Njj(O)  is  the  aerosol  number  density  at  ground  level  as 
determined  by  a  particle  counter  for  instance,  and  h  is 
the  altitude -in  kilometers. 

4.1.3    Ozone  Absorption 

We  now  consider  the  third  attenuating  component,  the 
ozone  absorption,  which  occurs  in  three  wavelength  bands i 
the  Hartley  band,  the  Huggins  band,  and  the  Chappius  band. 

The  Hartley  band  extends  from  200nm  to  300nm  but  is 
not  very  interesting  for  the  ground  level  observer  because 
the  optical  depth  is  so  large,  exceeding  100  at  255nm, 
that  no  photons  reach  the  ground  below  290nm. 

The  Huggins  band  is  of  greatest  interest  because  it 
is  the  dominant  attenuator  in  the  middle  ultraviolet. 
The  band  extends  from  300nm  to  360nm  decreasing  approx- 
imately exponentially  over  this  range.  If  the  structure  of 
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the  spectral  absorption  coefficient  is  not  important,  the 
following  simple  form  is  useful: 


where  A  is  in  nanometers  and  v;   is  the  equivalent  thick- 
ness  of  ozone  at  standard  temperature  and  pressure  expressed 
in  centimeters. 

The  Chappius  band  extends  from  440nm  to  850nm,  reach- 
ing a  maximum  at  600nm.  This  band  is  ordinarily  only  a 
small  part  of  the  total  optical  depth,  but  it  is  necessary 
to  explain  certain  twilight  phenomena.  In  Figure  4.2  we 
show  the  Huggins  and  Chappius  absorption  coefficients. 

The  ozone  attenuation  coefficient  is  given  by 


^oz^^)  =  Aoz<^)  ^oz  Noz^h)  (4.19) 


where-  the  Aoz  in  cm"  are  the  Vigroux  ozone  absorption 

coefficients  (32),,w   is  the  total  ozone  thickness  in  cm, 

and  N   in  km"  is  the  ozone  altitude  distribution.  N 

oz  oz 

is  also  represented  by  a  sum  of  generalized  distribution 
functions . 
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Ozone   cibsorption  coefficients   in  Huggins 
band    (300-360)    and  Chappius  band    (440-850) 


Figure  4.2 


42 


0.2139  exp{(h-23.0)/4.44} 
(  1  +  exp{(h-23.0)/4.44}  )' 


+  0.0096  e  h/5-78  (4.20) 

We  used  the  absorption  coefficients  of  Vigroux  for  -44 "C, 
but  Waltham  (33)  indicates  that  for  our  latitude  an 
average  ozone  layer  temperature  of  -55 "C  would  be  more 
accurate. 

The  optical  depths  of  all  three  major  attenuating 
components  are  shown  as  a  function  of  wavelength  in 
Figure  4.3,  the  standard  atmosphere  of  Elterman  being 
used  as  the  model. 

In  addition  to  the  attenuators  discussed  there  will  be 
absorption  by  oxygen,  water  vapor,  carbon  dioxide  and  other 
minor  gases.  None  of  these  absorb  light  in  our  wavelength 
range  except  H2O  which  has  bands,  which  we  avoid,  as  listed 
in  Table  4.1. 
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TABLE  4.1 


Water  vapor  absorption  bands. 


Designation    Wavelength  (nm)      Centroid  (niti) 


a 

700  -  740 

718 

SOOnm 

790  -  840 

810 

p,a,T 

930  -  980 

935 

«l> 

1095  -  1165 

1130 

* 

1319  -  1498 

1395 
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4.2    Standard  Bouguer-Langley  Analysis 

The  standard  technique  for  analyzing  solar  radiometer 
data  is  based  on  the  attenuation  law  attributed  variously 
to  Bouguer,  Beer,  and  Lambert. 


I  =  Iq  e""^  (4.21) 


where  the  total  optical  depth  is  the  sum  of  the  three  terms 
given  in  eqn.  4.1.  Unless  otherwise  stated  the  optical  depth 
is  for  ground  level  and  an  air  mass  of  one.  For  non-zero 
solar  zenith  angles  the  optical  depth  in  eqn.  4.21  is 
replaced  by  mx. 

The  turbidity  coefficient  used  here  is  that  of 
Volz  (34) .  The  turbidity  B  is  the  decadic  version  of  the 
aerosol  optical  depth  for  the  wavelength  500nm.  The  solar 
radiation  J  observed  by  a  Volz  sunphotometer  is  defined 
by  the  equation 

J*s  =  Jq  lO'^'^R  ^   "^oz  ■*•  ^^^  (4.22) 

where  s,  the  square  of  the  sun-Earth  distance  in  AU, 
corrects  for  variations  in  the  distance  and  Jq  is  the 
extraterrestrial  constant  for  the  instrument,  correspond- 
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ing  to  the  reading  the  instrument  would  give  at  the  top 
of  the  atmosphere. 

Sun photometers  are  frequently  calibrated  by  the 
Langley  method.  A  series  of  meter  readings  J  are  taken 
over  a  wide  range  of  air  mass  values  (m=l  to  m=6) .  If  the 
total  optical  depth  remains  constant  over  the  period  of 
observations,  and  this  requirement  is  the  weakness  of  the 
Langley  method,  then  a  graph  of  logJ  versus  m  will  be  a 
straight  line  and  extrapolate  to  Jq  at  m=0. 

Assuming  that  Jq  has  been  determined  and  that  the 
Rayleigh  and  ozone  base  10  optical  depths  at  500nm  are 
known  (  T^  =   0.0634  and  T   =  0.004)  equation  4.22  may  be 
solved  for  the  turbidity. 

B  =  i  log(  _  )  -  (  T  +  T   )  (4.23) 

J'S 

The  turbidity  was  originally  defined  for  500nm  only, 
but  the  definition  may  be  extended  to  other  wavelengths. 
The  turbidity  has  been  found  usually  to  vary  with  wave- 
length according  to 


B(A)  =  B(Aq)  (X/Aq)"°'  (4.24) 

The  parameter  a,  which  to  some  extent  characterizes  the 
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aerosol  size  distribution,  varies  from  0  to  4  with  a  world- 
wide average  of  1.3. 

4.3    Ultraviolet  Solar  Radiometry 

We  now  consider  certain  difficulties  encountered  in 
the  ultraviolet  and  not  in  the  visible.  The  equation  for 
the  signal  from  an  instrument  with  bandwidth  AX  centered 
at  \q  and  with  an  acceptance  angle  AQ   is 


J(X^)  =  C(A.)  /   /   R(X)  I{X,fl)  dX  dn        (4.25) 
"       "   An  AX 


where  J  is  the  observed  signal,  R(X)  is  the  relative 
spectral  response  of  the  instrument,  and  I{\,Q)    is  the  mono- 
chromatic intensity  at  wavelength  X  from  direction  fi. 
C(Xq)  is  the  calibration  constant  which  converts  the  signal 
to  absolute  irradiance. 


C(X„)  =  ^— y (4.26) 

/  R(X)  H^(X)  dX 
AX      ^ 


Here  H^  is  the  spectral  irradiance  of  a  standard  lamp. 

The  light  entering  a  sunphotometer  is  composed  of  the 
directly  transmitted  sunlight  and  the  scattered  sunlight. 
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1{X,Q)    =  I^{A,fi)  +  Ig(X,fi)  (4.27) 


Usually  the  instrument  acceptance  angle  is  one  or  two 
degrees  and  is  centered  so  as  to  contain  the  sun's  half 
degree  disc,  so  the  direction  integration  for  I,  is  not 
necessary.  The  scattered  intensity  is  variable,  but 
normally  it  is  so  much  less  intense  than  the  direct  compo- 
nent, by  three  or  four  orders  of  magnitude,  that  it  can  be 
neglected.  At  ultraviolet  wavelengths  scattering  is  much 
more  important.  But  aureole  measurements,  which  give  the 
ratio  of  scattered  intensity  to  direct  intensity,  show  that 
for  a  well  collimated  instrument  the  contribution  to  the 
total  signal  made  by  scattered  light  at  309nm  is  less  than 
1%  and  may  be  neglected. 

Inserting  the  solar  direct  intensity  at  ground  level 
into  eqn.  4.25  gives 

J(Xq)  =  C(X  )  /   R(A)  Hq(A)  e"^  '"i  ^i^^^  dA    (4.28) 
AA 

where  H„  is  the  extraterrestrial  solar  flux  and  the  m.  are 
the  individual  air  masses  for  the  three  attenuating  compo- 
nents. The  limit  m=0  gives  the  extraterrestrial  constant 
for  the  instrument. 
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Jo(Aq)  =  C(Xq)  /   R(A)  H  (A)  d\  (4.29) 

AX 


For  wavelengths  below  350nra  the  rapid  increase  and 
structure  in  the  ozone  spectral  absorption  coefficients 
lead  to  significant  change  in  the  spectral  quality  of 
attenuated  sunlight  with  air  mass,  producing  deviations 
from  straight  line  Bouguer-Langley  (BL)  plots.  Equation 
4.28  accurately  models  the  radiometer  to  allow  for  this 
effect.  Solar  irradiance  data  together  with  optical  depth 
data  were  input  at  O.lnm  intervals  into  eqn.  4.28.  A  family 
of  curves  is  generated  by  varying  Tj^^  and  plotting  the  log 
of  the  expected  instrument  reading  against  the  air  mass. 
Curves  generated  for  the  313nm  channel  are  shown  in  Figure 
4.4. 

A  form  of  Langley  extrapolation  to  find  J  is  possible 
using  such  a  family  of  curves  by  superimposing  data  for  a 
constant  turbidity  day  over  the  calculated  curves  and 
adjusting  Jq  to  find  the  best  fit.  This  method  gives  results 
as  good  as  calibration  with  a  standard  lamp. 

Because  of  deviations  from  straight  line  BL  plots  the 
simple  siinphotometry  equations  of  the  last  section  can  no 
longer  be  used;  yet  we  wish  to  be  able  to  calculate  the 
aerosol  optical  depth  from  air  mass  and  instrument  readings. 
Since  the  optical  depth  is  the  slope  of  the  straight  line 
of  a  BL  plot,  we  shall  consider  the  derivative  of  J(A«) 
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with  respect  to  m  (=  sec0)  and  designate  it  the  effective 
optical  depth. 


®        J{X«)     dm 


(4.30) 


or 


/  R{X)  Hq(A)  (  ^"^i  dm  J  e'^^i^^i 


dX 


T  (m)  = r-— (4.31) 

®  /  R(A)  lif.[X)    e"^^i"^i  dX 


The  dependence  of  t  (m)  on  m  is  shown  in  Figure  4 . 5  for 
four  values  of  the  aerosol  optical  depth.  The  degree  to 
which  the  curves  are  parallel  suggests  that  within  certain 
air  mass  limits  the  aerosol  optical  depth  is  independent 
of  air  mass. 


Tg(m)  =  TgQ(m)  +  T^  (4.32) 


T  Q (m)  is  T  (m)  calculated  for  zero  optical  depth  for 
the  aerosol. 

We  investigated  several  analytic  forms  for  t  q  and 
for  the  313nm  and  340nm  channels  foiand  the  following  form 
the  most  useful: 
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m  =  sec (6) 

Variation  of  total  effective  optical  depth 
with  aerosol  optical  depth  at  four  values. 


Figure  4 . 5 


53 


^0 
T^^dn)  =  (4.33) 

(  1  +  mt^)Y 

where  t^   is  t^^  at  m=0  and  t^  and  y   are  constants  for  a 
particular  channel  determined  using  a  least  squares  program 
to  fit  the  curve  generated  by  eqn.  4.31.  The  upper  curve 
in  Figure  4.6  shows  the  best  fit  for  the  313nm  channel, 
obtained  for  t^=0.6  86,  t^=0.293,  and  y=0.488. 

This  equation  is  integrable  when  inserted  into 
eqn.  4.30. 


InJ^CXg)  -lnJ(AQ)  =  mx^   ^ 


(1-Y)t3^ 
X  {  1  -  (l+mtj^)-'-"^}      (4.34) 

Solving  for  the  aerosol  optical  depth  gives  the  equivalent 
to  eqn.  4.23  which  we  were  seeking  for  use  in  the  ultra- 
violet. 


m 


mtj^d-Y) 
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For  the  330nm  channel  another  analytic  form  for  t  « 

eO 

was  found  useful. 


y+e  X 


This  is  the  lower  curve  in  Figure  4.6,  the  best  fit  being 
given  by  values  of  tQ=0.319,  t^=0.204,  and  y=6.35. 

When  we  investigated  the  error  involved  in  approxi- 
mating the  effective  optical  depth  as  the  siam  of  a  zero 
aerosol  term  containing  all  variation  with  m  plus  a  constant 
aerosol  optical  depth  term,  eqn.  4.32,  we  found  that 
the  error  was  small  for  the  313nm  channel:  about  3%  for 
an  air  mass  of  2.0.  But  we  found  that  the  error  for  the 
330nm  channel  was  almost  10%,  so  a  simple  correction  term 
for  the  residual  m  dependence  in  t^.  was  included  to  reduce 
the  error  to  less  than  1%. 


Tg(m)  =  TgQ(m)  +  Tj^(0.96  +  0.02m)     (4.37) 


This  form,  when  inserted  into  eqn.  4.30  and  integrated, 
gives  the  equation 
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■^M  = 


(0.96m-0.01m^) 


'    i„(^oiM.)  .  ^o"*^' 


I 


J{Aq)-s       y 


]  } 


(4.38) 


This  equation  was  used  for  analytic  calculation  of  aerosol 
optical  depths  for  the  330nin  channel. 

Clearly  the  equations  developed  here  are  for  our 
specific  instruments,  but  similar  equations  with  variations 
in  the  analytic  forms  as  required  will  be  applicable  for 
ultraviolet  radiometers  in  general. 

The  term  we  chose  to  study,  the  derivative  of  J(>^g) 

with  respect  to  m,  is  not  unique.  A  different  definition 

for  T  (m)  could  have  been  made  that  would  contain  all  m 
e 

variation  yet  not  require  re-integration  to  find  the  aerosol 
optical  depth.  Such  a  procedure  would  have  been  simpler  in 
some  ways,  but  it  might  have  made  the  analytic  fits  more 
difficult  to  achieve. 
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4 . 4    Instrumentation 

As  part  of  the  NOAA  (EPA)  national  turbidity  network 
(35)  we  have  made  routine  measurements  with  a  two  channel 
selenium  cell  sunphotometer.  The  channels  of  this  instru- 
ment are  centered  at  381nm  and  498nm  with  halfwidths  of 
llnm  and  13nm.  In  addition  we  have  a  six  channel  instrument 
kindly  provided  by  Dr.  F.  E.  Volz  which  has  channels  at 
313nm  and  341nm  with  halfwidths  of  12nm  and  13nm.  The 
two  visible  channels  overlap  the  first  instrument,  and 
we  have  not  made  use  of  the  two  infrared  channels,  vi^hich 
are  designed  to  observe  a  water  vapor  band. 

Selenium  cell  sunphotometers  are  widely  used  because 
of  their  simplicity,  low  cost,  and  portability.  They  do 
have  two  possible  drawbacks :  the  photocell  may  not  be 
completely  linear  throughout  its  range,  and  the  sensitivity 
and  spectral  response  of  the  instrument  may  change  with 
age.  One  can  correct  for  nonlinearity,  but  changes  in 
spectral  response  and  sensitivity  amount  to  a  change  in 
the  extraterrestrial  constant  of  the  instrument.  Such 
changes  necessitate  periodic  recalibration  of  the  instrument 
either  by  the  Langley  method  or  by  comparison  with  an 
infrequently  used  standard  instrument.  We  have  observed 
some  drift  in  our  49  8nm  and  381nm  channels  over  the  past 
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year,  while  the  313nm  and  341nin  channels  have  been  fairly 
stable. 

In  addition  we  have  constructed  a  research  quality 
solar  radiometer  for  use  in  conjunction  with  our  field 
instruments.  The  instrument  was  first  used  with  an  RCA  7102 
photomultiplier  tube  operated  at  680  volts.  This  is  a  low 
gain  photomultiplier  tube  which  has  as  S-1  photocathode 
of  maximum  sensitivity  at  850nm  but  still  useful  at  330nm. 
The  photometer  can  usefully  span  three  decades  in  intensity 
and  is  linear  to  within  3%  over  this  range  (36)  .  It  is 
used  in  conjunction  with  interference  filters  of  lOnm  half- 
width  centered  at  wavelengths  of  330nm,  400nm,  500nm,  600nm, 
and  700nm.  The  field  of  view  of  the  instrument  is  limited 
by  a  two  aperture  collimator  with  an  acceptance  angle  of 
68*  of  arc,  sufficient  to  contain  the  32'  solar  disc  while 
admitting  little  scattered  light. 

Recently  this  instrument  was  modified,  replacing  the 
photomultiplier  tube  with  an  RCA  1P28  which  has  an  S-5 
photocathode  of  maximum  sensitivity  at  350nm  and  useful  at 
700nm.  Due  to  damage,  the  400nm  and  500nm  interference 
filters  were  replaced  with  420nm  and  505nm  absorption 
filters  with  bandwidths  of  50nm  and  30nm.  Such  large  band- 
widths  would  normally  be  avoided,  but  use  of  eqn.  4.2  8 
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Multichannel  photomultiplier  solar  radiometer 
mounted  alongside  camera  emd  occulting  disc 
used  to  photograph  the  solar  aureole. 


Figure  4,7 
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to  model  these  channels  allowed  us  to  use  them  with  no 
problems.  In  addition  the  five  filters  were  installed  in  a 
motor  driven  filter  wheel,  allowing  rapid  observations  to 
be  made  and  use  of  a  chart  recorder  if  desired.  The  photo- 
multiplier  solar  radiometer  is  shown  in  Figure  4.7. 

4 . 5    Results 

We  made  several  series  of  observations  in  1974,  and 
of  these  three  days  exhibited  sufficiently  constant  aerosol 
optical  depth  for  the  BL  method  to  be  useful:  February  4, 
February  10,  and  November  22.  The  February  observations 
were  made  in  developing  our  analytic  technique  for  calcu- 
lating ultraviolet  aerosol  optical  depths,  while  the 
November  data  were  taken  in  conjunction  with  the  solar 
aureole  study  discussed  in  Chapter  5. 

These  three  days  were  all  characterized  by  very  low 
optical  depths  as  an  anticyclonic  mass  of  cold  Cemadian 
air  moved  into  Florida.  The  more  usual  atmospheric  condition 
is  a  low  to  medium  dust  level  with  water  droplet  haze, 
because  of  the  warm  hvunid  air.  This  leads  to  a  very  varia- 
ble aerosol  optical  depth.  A  common  pattern  might  be  a  low 
morning  optical  depth  of  0.1  increasing  to  0.2  by  noon, 
then  decreasing  somewhat  in  the  afternoon. 
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In  Figures  4.9  through  4.12  we  show  photometer  readings 
as  a  function  of  air  mass  for  various  of  these  days  along 
with  calculated  BL  plots  for  various  aerosol  optical  depths 
generated  by  our  computer  modeling  of  the  instruments. 

In  Figure  4.9  we  show  readings  from  the  700nm  and 
600nm  channels  for  February  4  and  November  22,  1974.  It 
shows  that  the  optical  depth  for  November  22  was  about  half 
that  of  February  4,   and  that  while  Tj.  at  600nm  was  as 
expected  slightly  larger  than  at  700nm  for  the  February 
data,  for  the  November  data  it  was  slightly  larger  at 
700nm  than  at  600nm. 

Data  for  the  505nm  and  420nm  channels  shown  in  Figure 
4.10  are  given  for  November  only,  since  the  new  filters 
were  only  recently  installed.  These  readings  followed 
straight  lines  very  closely,  ideal  for  Langley  extrapolation 
to  Jfj.  From  420nm  to  600nm  t„  decreases  with  wavelength 
approximately  as  a=1.5. 

Plotted  along  with  the  505nm  data  is  a  set  of  data 
taken  with  the  Volz  six  channel  instrument,  500nm  channel, 
on  April  27,  1975  in  Cedar  Key,  Florida,  an  island  off  the 
west  coast  of  Florida.  An  optical  depth  of  0.32  is  indicated 
for  Cedar  Key  for  this  day.  The  point  plotted  at  m=1.05 
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Data  plotted  with  computer  generated  Bouguer- 
Langley  plots  for  x  =  0.0,  0.046,  0.115,  0.23, 
0.46,  and  0.69. 


Figure  4.10 
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is  a  measurement  made  in  Gainesville  on  the  same  day. 
T^   for  Gainesville  was  0.60,  almost  twice  that  at  the 
coast.  This  is  due  partly  to  lower  dust  levels  over  the 
ocean,  but  more  important  is  the  solar  heating  of  land 
leading  to  thermals  and  a  high  rate  of  cloud  (and  haze) 
formation. 

Figure  4.11  shows  February  4  data  for  the  330nm 
channel.  On  this  particular  day  the  data  follow;  the 
nonlinear  calculated  BL  curves  quite  well. 

In  Figure  4.12  data  for  the  313nm  channel  for  two 
days,  February  10  and  November  22,  are  shown.  The  point 
at  m=4.7  which  falls  above  the  t  =0  line  is  probably  due 
to  inaccurate  zeroing  of  the  instrument,  since  an  error 
equal  to  the  smallest  scale  division  of  the  instrument  is 
indicated. 

On  November  22  the  aerosol  optical  depths  for  the 
700nm,  600nm,  505nm,  420nm,  and  313nm  channels  respective- 
ly were  0.030,  0.018,  0.023,  0.035,  and  0.12.  Except  for 
the  700nm  channel,  this  follows  the  expected  decrease  of 
optical  depth  with  wavelength  ,  with  alpha,  the  wavelength 
exponent  defined  by  eqn.  4.24  in  the  range  1.2  to  1.6. 
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CHAPTER  5 


THE  SOLAR  AUREOLE 


The  light  that  constitutes  the  daytime  sky  brightness 
is  composed  partially  of  light  that  has  been  Rayleigh 
scattered  and  partially  of  light  that  has  been  Mie  scat- 
tered. If  one  is  interested  in  determining  characteristics 
of  the  atmospheric  aerosol  by  studying  light  that  has  been 
Mie  scattered,  it  is  logical  to  concentrate  on  that  region 
of  the  sky  in  which  Mie  scattering  is  maximum:  the  solar 
aureole.  The  aureole  is  the  circular  area  of  whitish  light 
around  the  sun,  ten  or  so  degrees  in  extent,  arising  from 
the  sharply  forward  peaked  scattering  pattern  from  aerosol 
particles. 

5.1    Aerosol  Models 

The  variability  of  the  natural  aerosol  was  discussed 
in  Chapter  2,  and  it  may  be  concluded  that  the  size  distri- 
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bution  is  the  important  variable  in  solar  aureole  study. 
In  its  most  general  form  the  size  distribution  includes 
variation  with  altitude,  but  using  a  size-altitude  distri- 
bution adds  greatly  to  the  complexity  and  cost  of  an 
atmospheric  scattering  calculation  while  yielding  little 
increase  in  accuracy. 

The  most  prominent  feature  of  the  natural  aerosol  size 
distribution  is  the  power  law  decrease  in  number  density 
v/ith  increasing  particle  radius  first  described  by  Junge 
(10) .  Here  the  term  size  distribution  refers  to  the  radius- 
number  distribution  n(r) :  the  number  of  particles  per  cubic 
centimeter  within  a  unit  radius  range  dr  at  x,   where  r  is 
the  particle  radius  in  microns.  Another  useful  function  is 
the  cumulative  oversize  distribution  N(r) :  the  total  number 
of  particles  per  cubic  centimeter  that  have  radii  greater 
than  r.  Thus  N(0)  will  be  the  total  number  density  of 
particles.  The  size  distribution  is  related  to  the  cumula- 
tive oversize  distribution  by 


N(r)  =  /^  n(r)  dr  (5.1) 


or  alternatively 
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n(r)  =  -dN(r)/dr  (5.2) 

Junge  prefers  to  use  a  logarithmic  scale  because  of  the 
wide  range  of  particle  sizes  and  concentrations.  His  log 
radius-number  distribution  is  defined  as 

Nj^(r)  =  -dN(r)/d(log  r)  (5.3) 

He  found  that  continental  aerosols  are  described  by 


dN(r)  =  c  r""^  d(log  r)  (5.4) 


or 


n(r)  =  0.434  c  r"^^"*"^^  (5.5) 


This  distribution  is  frequently  used  directly  in  calcu- 
lations with  an  externally  imposed  lower  limit  (  about  0.01 
microns)  to  avoid  r  =  0  problems.  The  constant  c  depends 
on  the  number  of  particles  per  cubic  centimeter  and  the 
exponent  v  is  found  to  range  between  2.5  and  4.0.  The 
majority  of  aerosol  measurements  indicate  that  v  =  3 
gives  a  good  fit  to  the  natural  aerosol  distribution. 
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TABLE   5.1 


Parameters    for  Deirinendjian' s   size  distributions, 


HAZE  MODEL 

PARAMETER 

M 

L 

H 

a 

5. 333x10** 

4.976x106 

4.000x10^ 

b 

8.9443 

15.1186 

20.0 

a 

1.0 

2.0 

2.0 

Y 

0.5 

0.5 

1.0 

max 

O.OSp 

0.07y 

O.lOy 
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A  very  widely  used  functional  form  for  the  size 
distribution  is  that  of  Deirmendjian  (37) 


n(r)  =  a  r"  e"  ^^  ^^^  (5.6) 


which  vanishes  at  r  =  0  and  which  he  called  a  modified 
gamma  function,  since  it  reduces  to  the  gamma  function  when 
Y  =  1.0.  The  parameters  for  three  models  are  listed  in 
Table  5.1.  Haze  M  was  introduced  to  reproduce  a  typical 
marine  aerosol,  while  L  represents  a  continental  aerosol. 
Haze  H  represents  the  stratospheric  aerosol  submicron 
particles. 

A  number  of  models  were  tried  by  Green,  Deepak,  and 
Lipofsky  (36)  ,  the  most  useful  being  their  mathematical 
spline  model 


exp(f]''/  [{bXexpCfj'^-D+l}^  r^" 


n(r)  =  c  exp(|]V  ^{b(exp(f)  ^-11  •  -  ^  -^+3 

X  Q(r) 


(5.7) 


where  a,  b,  c,  and  v  are  adjustable  constants  and  Q(r) 
is  the  efficiency  factor. 

In  our  analysis  of  light  scattering,  we  use  a  regular- 
ized power  law  that  is  expressed  simply  in  teinns  of  the  over- 
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size  distribution 


N(r)  =  N(0)/(  1  +  (|)^)  (5.8) 


which  corresponds  to  the  size  distribution 

n(r)  =  N(0)  fv  ^-= (5.9) 

This  size  distribution  also  vanishes  at  r  =  0 ,  and  for 
r  >>  a 

n(r)  ^  N(0)  ^^   ^   2v  (5-10) 


^  N(0)  ^^  r-<^^l)  (5.11) 

a 

which  is  identical  with  the  Junge  distribution.  This  para- 
meter V  thus  corresponds  to  that  of  Junge,  and  a  is  approx- 
imately the  radius  at  which  the  number  frequency  is  maximum. 
Junge  (38)  was  one  of  the  first  to  obtain  a  "complete" 
size  distribution  curve  in  1955  using  an  impactor  to  deter- 
mine the  distribution  of  particles  larger  than  0.05  microns 
and  an  ion  counter  to  size  the  particles  smaller  than  0.05 
microns.  Figure  5.1  shows  our  regularized  power  law 
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distribution  for  three  different  values  of  v,  and  in 
Figure  5.2  we  compare  the  data  points  found  by  Junge  with 
Deirmendjian's  standard  Haze  L  distribution  and  our  regu- 
larized power  law  distribution  adjusted  to  v  =  2.8  and 
a  =  0.03.  Pasceri  and  Friedlander  (39)  support  this  type 
of  distribution,  peaking  between  0.01  micron  and  0.10  micron 
and  then  decreasing  as  a  power  law,  with  more  recent 
measurements  made  in  1965  in  Baltimore  using  a  cascade 
impactor  and  a  rotating  disc  sampler. 

While  the  simple  regularized  power  law  is  adequate 
to  describe  a  normal  size  distribution,  an  additional  term 
can  be  added  to  describe  a  bimodal  distribution. 

n(r)  =  N^  \   ^ +  N.  -J   ^-— (5.12) 

'a-  a.(|)^)a     ^a^   (l.(|J^^)^ 

where  the  second  term  may  be  negative  if  required. 

We  attempt  to  infer  the  size  distribution  based  on 
the  scattering  pattern  of  the  solar  aureole  and  the  aerosol 
optical  depth  as  measured  by  the  solar  radiometer,  yet 
this  represents  only  a  limited  amo\int  of  information.  Our 
technique  is  to  constrain  the  aerosol  size  distribution 
models  to  those  that  are  physically  reasonable  in  light  of 
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the  cumulation  of  prior  investigations.  Other  investigators 
handle  the  problem  of  inversion  to  find  the  size  distri- 
bution without  assuming  a  size  distribution  model.  The 
Backus-Gilbert  inversion  technique  (40)  for  instance  may 
be  capable  of  establishing  the  size  distribution  with  no 
a  priori  assumptions.  But  it  has  generally  been  found  that 
noisy  data  present  difficulties  to  strict  numerical  inver- 
sion methods.  Ours  is  a  parametric  model-comparative  test 
with  data  approach  which  has  been  used  successfully  in  a 
number  of  fields,  particularly  in  nuclear  and  particle 
physics  where  scattering  measurements  have  reached  a  high 
level  of  sophistication. 

5.2    Mie  Phase  Fxonctions 

The  phase  function  relates  the  particle  size  distri- 
bution  (in  the  case  of  Mie  scattering)  to  the  intensity 
of  light  scattered  as  a  function  of  angle.  In  the  case  of 
Rayleigh  scattering  the  phase  function  is  analytic 


^R^^^  =  ri-?  (  1  +  °°s'e  )  (5.13) 


and  satisfies  the  normalization  condition 

^rC  ^^^^  ^^^^   ^^  <^*  =  1  (5.14) 
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Mie  scattering  intensities  are  calculated  in  terms 
of  Bessel  fimctions  and  are  not  simply  analytic,  so  the 
phase  function  which  is  the  result  of  averaging  Mie  scat- 
tering over  an  arbitrary  size  distribution  certainly  is 
not  analytic.  The  basic  definition  for  a  phase  function  is 
that  it  is  the  ratio 

a„{e) 

P(e)  =  2 (5.15) 

/  o„(e)  do) 

of  the  differential  scattering  cross  section  to  the  total 
cross  section.  In  terms  of  the  complex  amplitudes  for  the 
scattered  radiation  and  the  scattering  efficiency  factor 
of  Van  de  Hulst 

^  2Tr  ^      2  TT  r^  Q 

s 

This  phase  function  for  a  single  particle  is  then  averaged 
over  the  size  distribution 


If  n(r)  is  our  regularized  power  law  distribution,  the 
phase  functions  at  different  wavelengths  are  related  by 
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Pj^(e,v,a,A)  =  Pj^(e,v,ka,kX)         (5.18) 

so  that  once  the  phase  function  is  calculated  at  one  wave- 
length for  an  array  of  v  and  a  values,  it  is  directly 
extensible  to  a  range  of  wavelengths. 

Since  we  infer  the  size  distribution  by  varying  the 
model  parameters  v  and  a  to  find  the  best  correspondence 
with  experiment,  it  would  be  more  efficient  as  a  matter  of 
practical  search  procedure  if  we  could  characterize  the 
phase  function  analytically  as  a  function  of  angle.  Then 
we  could  search  directly  on  variations  of  the  phase  function 
and  save  the  time  and  expense  of  calculating  new  phase 
functions  for  every  variation  of  the  size  distribution. 

Henyey  and  Greenstein  (41)  proposed  the  scattering 
function 

P(e)  =  x"^/V4Tr  (5.19) 

where 

X  =  1  +  g2  -  2g  cose  (5.20) 

Dr.  A.E.S.  Green  has  proposed  a  generalized  Henyey- 
Greenstein  function 
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g'T  r       1  1       1-1 


27r  x^+^/2  '^      |l-gr       (1+g   ) 


1-1 I—] 

^       ll-qT        (1+q    )■'    J 


(5.21) 


which  contains  the  Henyey-Greenstein  function  as  a  special 
case  when  t  =  1.  The  interesting  domain  of  sharply  forward 
peaked  phase  functions  is  characterized  by  values  of  g 
near  1. 

Recognizing  the  similarity  of  sharply  peaked  light 
scattering  functions  to  the  strong  forward  peaking  for 
allowed  transitions  in  inelastic  scattering  of  high  energy 
electrons  by  atoms.  Green  developed  the  Born-Bethe  approx- 
imation analogy.  In  this  approximation  the  differential 
cross  section  is  given  by 


do  _  4  k'  f (x)  fc   o,v 

dfl    "    w    k    ~ir~  ^^"^^' 


where  k  is  the  incident  electron  energy  and  k'  the  final. 
If  W  is  the  excitation  energy,  x  is  the  square  of  the 
momentum  transferred  and  is  given  by 


X 


=  k^  +  k'^  -  2kk'  cose  (5.23) 


One  of  the  powerful  features  of  the  Born-Bethe  approx- 
imation is  that  if  one  knows  the  angular  distribution  at 
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one  energy,  it  can  be  found  at  any  other  energy,  closely 
akin  to  what  we  are  trying  to  do  with  regard  to  wavelength. 
Recalling  eqn.  5.20  we  see  that  the  analogues  k=l  and 

k'=g  places  the  two  x's  in  correspondence.  Introducing  the 

2 
parameter  w  =  1-g  ,  it  follows  that 


X  =  (  1  -  /  1-w) ^  +  2  /  1-w  (  1-cose  )      (5.24) 

a  complex  relationship  which  can  be  simplified  by  introduc- 
ing the  parameter  v 


V  =  2(  1  -  /~r=w   )  (5.25) 


in  terms  of  which  the  phase  function  may  be  written  as 


P(e)  =  I  (5.26) 


where 

2 


x=j+   (2-v)(l-  cose  )        (5.27) 


Better  fits  to  the  actual  phase  functions  were  obtained 
for  the  scattering  range  0"  to  20°  by  adding  a  constant 

P(e)  =  I  +  B  (5.28) 
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TABLE    5.2 


Phase   ftmction  parameters  v,  A,   and  B  as   a  function 
of  V   and  a. 


2.0 

2.5 

V 

3.0 

3.5 

4.0 

0.01 

.0203 
.0203 
.2440 

.0286 

.00968 

.655 

.0437 

.00444 

.741 

.0766 

.00283 

.627 

.1480 

.00304 

.477 

0.03 

.0203 
.0211 
.299 

.0286 
.0102 
.664 

.0438 

.00475 

.768 

.0768 

.00310 

.664 

.1490 

.00347 

.520 

0.05 

.0202 
.0223 
.201 

.0286 
.0112 
.670 

.0440 
.0054 
.808 

.0772 

.00368 

.722 

.1480 

.00430 

.584 

a 

0.10 

.0200 
.0260 
.0910 

.0286 
.0149 
.622 

.0445 

.00835 

.870 

.0782 

.00643 

.849 

.1450 

.00791 

.731 

0.30 

.0187 
.0366 
-.332 

.0269 
.0324 
-.0105 

.0420 
.0292 
.243 

.0697 
.0314 
.291 

.1100 
.0392 
.186 

0.50 

.0172 
.0396 
-.516 

.0239 
.0407 
-.394 

.0353 
.0416 
-.287 

.0537 
.0463 
-.297 

.0783 
.0552 
-.415 
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In  Table  5.2  we  give  the  best  fit  v.  A,  and  B  para- 
meters for  an  array  of  v  and  a  size  distribution  phase 
functions.  For  an  arbitrary  v  and  a  combination,  interpo- 
lation from  the  table  gives  the  three  parameters  needed  to 
generate  the  analytic  phase  function.  This  represents  a 
substantial  saving  in  computational  time,  since  direct 
generation  of  the  phase  function  involves  integrating 
the  Mie  scattering  coefficients  over  the  size  distribution. 

5.3    The  Single  Scattering  Calculation 

A  theoretical  calculation  of  sky  brightness  assuming 
that  light  reaching  the  ground  has  scattered  one  time 
only  produces  quite  good  agreement  with  observation  under 
conditions  of  small  aerosol  optical  depth  and  aureole 
scattering  angles.  In  previous  sections  we  have  developed 
the  basis  needed  for  the  single  scattering  calculation: 
the  altitude  distribution  models  for  air,  particulates, 
and  ozone,  the  absorptive  behavior  for  ozone  and  the  scat- 
tering behavior  for  aerosols  and  air,  and  the  models  to 
be  used  for  the  aerosol  size  distribution.  We  can  now  cal- 
culate expected  aureole  scattering  using  these  in  a  real- 
istic model  of  the  atmosphere. 

The  intensity  I  of  monochromatic  light  observed  in 
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the  direction  16, (p),    zenith  angle  and  azimuth  angle  of 
line  of  sight,  when  the  sun  is  at  zenith  angle  6  is 
given  by 


xIq  T^   T^    ds  (5.29) 


where  i|;(0,(j))  is  the  scattering  angle,  h  is  the  altitude 
during  integration  over  path  s,  and  Iq  is  the  solar 
intensity.  T,  and  T-  are  the  transmission  factors  from  the 
scattering  volume  to  the  ground  and  from  the  top  of  the 
atmosphere  to  the  scattering  volume. 


T2(h)  =  exp(-/"  {3j^{s)+3j,{s)+6Q^(s)}ds  )    (5.30) 


T^ih)    =   exp(-/^  {3j^(l)+3j^(l)+6^^(l)}dl  )    (5.31) 

Here  s  is  the  path  followed  by  the  sunlight  in  reaching  the 
scattering  volume,  and  1  is  the  path  (  of  length  L) 
followed  in  reaching  the  ground.  These  paths  should  include 
the  bending  of  light  due  to  atmospheric  refraction  for 
solar  zenith  angles  near  ninety  degrees,  but  we  did  not 
include  it  in  our  program  since  all  calculations  were  for 
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zenith  angles  less  than  ninety  degrees. 

The  altitude  distributions  used  in  our  model  atmos- 
phere were  the  analytic  functional  forms  described  in 
Chapter  4.  The  relative  altitude  distributions  are  assumed 
invariant;  changes  in  aerosol  optical  depth  and  ozone 
thickness  are  made  by  varying  the  altitude  distribution 
uniformly. 

5.3.1    Geometry 

As  noted  earlier  the  atmosphere  may  be  regarded  as 

plane  parallel  for  zenith  angles  of  up  to  75">.  Thus,  in  the 

small  angle  form  of  our  program  the  geometry  of  Figure 

5.3a  v/as  used.  If  the  sun  is  at  zenith  angle  6  and  the 

s 

point  P  in  the  sky  is  at  zenith  angle  6  and  relative 
azimuth  angle (})  to  the  sun,  sunlight  must  scatter  through 
an  angle  ^   to  be  observed  in  the  direction  of  P. 

cosT|j  =  cos  e  cos  e  +  sin  9  sin  9  cos  (J)   (5.32) 

s  s 

The  geometry  for  near-horizon  light  scattering  is 
diagramed  in  Figure  5.4.  We  simplify  the  problem  for  now 
by  restricting  the  calculation  to  points  along  a  vertical 
cut  through  the  sun  (lie.  along  a  line  of  constant  azimuth) . 
The  computer  integrates  along  the  path  L  from  groiuid  level 
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a)  Scattering  geometry,  plane  parallel 
atmosphere . 


b)  Angle  relationships, 


Figure  5 . 3 
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Spherical  atmosphere  scattering  geometry. 


Figure  5 . 4 
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to  altitude  A,  where  A  is  expressed  in  terms  of  the  vertical 
height  h  by 


A  =  (R^  +  2Rh  +  h^sec^e)^/^   -   R         (5.33) 


where  R  is  the  radius  of  the  Earth.  For  each  point  along 
L  sunlight  inust  travel  along  a  path  s,  for  which  the 
altitude  A'  is  given  by 


A'  =  [(R+A)^  +  (h'-h)2sec^eg  +  2(R+A)(h'-h) 


xsecS  cos{  0  -arcos  (pxr-)  } 


1/2 


-  R      (5.34) 


The  integration  is  carried  out  in  eighty  equal  increments 
in  h  to  a  maximum  altitude  A  of  eighty  kilometers.  We 
found  this  integration  step  size  to  give  excellent  results 
when  used  with  Simpson's  rule,  since  doubling  the  number 
of  integration  points  caused  only  a  0.01%  variation. 

Although  we  did  not  choose  to  look  at  non-zero  ij) 
angles,  the  generalization  goes  as  follows. 

The  integration  along  L  is  unchanged,  but  along  s 

A'  =  {(R+h)2  +  2(R+h)A  +  A^sec^Y)"''^^      (5.35) 

-R 
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where 

^°^Y  =  i^^ose^  -f  {l-(^)^}^/23ine  cos,f     (5.36) 


and 


^2  =  exp(-/"  B(A')secY  dh*  ]  (5.37) 

A 


The  scattering  angle  ip  will  still  be  given  correctly  by 
eqn.  5.32. 

Our  objective  in  performing  this  calculation  is  to 
compare  calculated  aureole  intensities  with  aureole  inten- 
sities obtained  by  our  photographic  technique,  which  is 
described  in  detail  in  Chapter  6.  An  example  of  an  aureole 
photograph  is  shown  in  Figure  5.5.  A  key  feature  of  the 
photographic  method  is  that  the  solar  disc  (attenuated  by 
a  neutral  density  4.0  filter)  is  recorded  along  with  the 
aureole,  and  it  is  more  accurate  to  obtain  the  ratio 
of  aureole  to  solar  disc  intensities  than  to  attempt  to 
obtain  absolute  intensities  directly  from  the  photograph. 
Once  this  ratio  is  accurately  known,  reference  to  the 
direct  solar  intensity  at  ground  level  measured  by  the 
solar  radiometer  allows  the  absolute  aureole  intensities 
to  be  inferred.  The  intensity  of  light  for  the  sun's  disc 
is  calculated  by 
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I«(X,e^)  =  -i-  e^'^oSecSg  (5.38) 

s 
where  Tg  is  the  total  optical  depth  at  ground  level  and 
An  is  the  solid  angle  sxibtended  by  the  sun. 


5.3.2 Wavelength  Integration 

The  aureole  data  obtained  experimentally  are  not 
monochromatic,  since  an  interference  filter  is  used  at 
309nm  and  narrow  band  absorption  filters  are  used  at  430nm, 
540nm,  and  640nm  (Wratten  filters  #98,  #99,  and  #92). 
Initially  the  complete  single  scattering  calculation  was 
performed  at  lOnra  intervals  (at  Inra  intervals  for  the 
309nm  data)  and  integrated  to  obtain  the  expected  aureole 
pattern. 

1(0  ,6)    /  Tp(A)Sp(x)i^(e3,e)  dA 

^     -   ^^  (5.39) 


^s^^s^      /   Tj,(A)Sp(X)Ig(0g,8)  dX 

A  A 

where  T  is  the  filter  transmission  and  S_  is  the  relative 
film  sensitivity.  But  aureole  scattering  does  not  change 
rapidly  with  wavelength,  so  we  now  calculate  the  scattering 
at  three  wavelengths,  the  filter  minimum,  peak,  and  max- 
imum wavelengths,  and  interpolate  the  scattered  intensity 
at  each  scattering  angle  to  wavelengths  in  between  for 
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integration.  The  extraterrestrial  solar  intensity,  which 
is  irregular  with  wavelength,  is  removed  before  the 
quadratic  interpolation  is  performed  and  then  re-inserted 
at  each  wavelength  after  interpolation.  This  saves 
substantial  computation  while  causing  less  than  1%  change. 

5.3.3    Ground  Albedo  Correction 

Since  light  traveling  at  a  zenith  angle  close  to 
ninety  degrees  is  close  to  the  ground  over  a  large  part  of 
its  path,  we  felt  that  including  light  reflected  from  the 
ground  would  increase  the  accuracy  of  the  near-sunset 
aureole  calculation.  Bauer  and  Button  (42)  measured  ground 
albedo  in  Wisconsin  from  an  airplane  for  different  seasons 
and  terrain  and  found  that  wooded  terrain  in  summer  had 
an  average  albedo  of  0.15  to  0.20.  The  flat,  sandy  soil 
of  Florida  might  be  expected  to  have  a  slightly  higher 
albedo,  so  we  used  a  value  of  0.25  in  our  calculation  to 
correspond  with  that  used  by  Coulson,  Dave,  and  Sekera  (43) 
in  their  Rayleigh  scattering  calculation.  We  assume  that  the 
ground  is  a  Lambertian  reflector,  an  assumption  supported 
by<  but,  as  they  are  careful  to  point  out,  not  required  by 
the  measurements  of  Bauer  and  Dutton. 

For  the  purpose  of  calculating  the  albedo  contri- 
bution, the  ground  beneath  the  scattering  volume  along 
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Scattering  geometry  for  grotind  reflection 
contribution. 


Figure  5.6 
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path  L  is  assumed  to  be  plane  parallel.  The  scattering 
geometry  is  shown  in  Figure  5.6.  The  area  element  da 
reflects  direct  sunlight  to  the  scattering  element  along 


r. 


cose  A^   da 
dl  =  IgCe^)  cose^  ^—^ '^g^^q^       (5.40) 

Trr 


where  I^  is  the  intensity  of  sunlight  at  area  element  da, 
Aq  is  the  ground  albedo,  and  T   is  the  transmission  along 
r,   calculated  as  in  eqn.  5.31.  The  contribution  to  the 
aureole  from  ground  reflection  is  found  by  integrating 
over  the  entire  plane. 


di^g(0g,e,h)  =  ig(eg,x)cosegTj^(h)sec6  rjj^i^^ih) 

^P^ir)    +  3^(h)P^(;p')}Tg(0^)sinegdegd4,g    (5.41) 


where  ^'    is  calculated  as  in  eqn.  5.32.  This  contribution 
is  added  to  the  integrand  in  eqn.  5.29. 

This  procedure  was  used  in  our  calculations  with  an 
Aq  of  0.25  except  for  the  309nm  calculation,  for  which  an 
albedo  of  0.0  was  used.  The  average  albedo  for  the  300nm 
region  is  0.05  approximately,  and  Shettle  (44)  found  that 
using  0.0  instead  of  0.05  caused  only  a  1%  change  or  less 
in  his  multistream  multiple  scattering  calculation.  We 
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found  the  contribution  from  ground  reflection  to  be  more 
than  a  couple  of  per  cent  only  for  very  clear  air;  it  was 
not  as  important  at  large  zenith  angles  as  expected.  For 
a  pure  Rayleigh  atmosphere  ground  reflection  can  contri- 
bute as  much  as  25%  in  increased  light,  but  for  the  compu- 
tations for  the  two  days  studied  here,  only  abov.t  a  3% 
increase  was  observed. 

5.3.4 Rayleigh  Multiple  Scattering  Correction 

The  major  process  neglected  in  any  single  scattering 
calculation  is,  of  course,  multiple  scattering.  Multiple 
scattering  is  important  at  short  wavelengths,  large  optical 
depths,  large  scattering  angles,  and  large  zenith  angles. 
Accurately  representing  multiple  scattering  requires  the 
use  of  a  multistream  calculation  such  as  that  of  Mudgett 
and  Richards  (45)  or  a  Monte  Carlo  calculation  such  as 
that  used  by  Plass  and  Kattawar  (7) .  Since  our  objective 
is  to  model  atmospheric  scattering  as  closely  as  possible 
without  resorting  to  a  full  multiple  scattering  formalism, 
we  add  an  empirical  term  corresponding  to  the  multiple 
scattering  from  a  pure  Rayleigh  component  only.  CouIsoq, 
Dave,  and  Sekera  (43)  have  published  tables  giving  the 
scattered  radiation  including  all  multiple  scattering  based 
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on  the  exact  solution  of  Chandrasekhar  (5).  Their  tables 
give  the  scattered  radiation  for  a  plane  parallel  Rayleigh 
atmosphere  for  various  solar  zenith  angles,  scattering 
angles,  albedos,  and  wavelengths. 

For  comparison  we  reproduced  the  tables  for  single 
scattering  only  in  order  to  see  the  effects  of  multiple 
scattering  on  the  Rayleigh  component.  In  Figure  5.7  we 
compare  the  multiple  scattering  results  and  the  single 
scattering  results  for  two  different  sun  angles  and  two 
different  wavelengths.  One  notes  that  the  multiple  scat- 
tering and  single  scattering  sky  brightness  patterns  are 
veiry  similar,  the  difference  being  primarily  one  of  magni- 
tude. To  a  good  approximation,  then,  the  results  of  the 
Coulson,  Dave,  and  Sekera  (CDS)  calculation  can  be  matched 
by  multiplying  the  single  scattering  result  by  a  constant 
correction  factor.  If  the  correction  factor  is  calculated 
for  a  scattering  angle  equal  to  iero  and  applied  only  to 
the  solar  aureole,  the  fifteen  degrees  or  so  in  the  vicinity 
of  the  sun,  then  the  Rayleigh  component  should  be  corrected 
fairly  accurately  for  multiple  scattering. 

Table  5.3  gives  the  CDS  correction  factors  for  an 
array  of  wavelengths  and  solar  zenith  angles.  The  variation 
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Rayleigh  multiple  scattering  as  calculated  by 
Coulson,  Dave,  and  Sekera  compared  with  single 
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Figure  5 . 7 
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TABLE  5.3 


Ratio  of  Rayleigh  intensity  according  to  Coulson, 
Dave,  and  Sekera  to  single  scattering  intensity 
for  the  "near  sun"  area  of  sky  as  a  fxonction  of 
wavelength  and  solar  zenith  angle. 


X 

0 

23.1 

36.9 

s 
53.1 

66.4 

78.5 

84.3 

312.0 

2.225 

2.277 

2.353 

2.556 

3.047 

6.778 

- 

371.5 

1.623 

1.637 

1.660 

1.712 

1.786 

2.010 

3.653 

436.5 

1.323 

1.331 

1.342 

1.369 

1.401 

1.453 

1.607 

495.0 

1.216 

1.218 

1.224 

1.240 

1.258 

1.275 

1.306 

546.0 

1.135 

1.144 

1.150 

1.162 

1.176 

1.186 

1.201 

644.0 

1.058 

1.061 

1.062 

1.071 

1.078 

1.082 

1.095 

809.0 

1.021 

1.011 

1.011 

1.022 

1.029 

1.024 

1.035 
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of  the  CDS  factor  with  wavelength  is  smooth  and  is  des- 
cribed by  the  function 

F(eg,X)  =  B  +  (C-B)xf5i6^j  (5.42) 

where  a,  B,  and  C  are  constants  for  a  given  solar  zenith 
angle  determined  using  a  nonlinear  least  squares  routine. 
The  dependence  of  F  (or  rather  of  a,  B,  and  C)  on  solar 
zenith  angle  is  not  so  regular,  so  the  three  constants 
for  an  array  of  angles  were  included  in  numerical  form. 

The  importance  of  the  CDS  factor  depends  on  the 
relative  optical  depth  of  the  aerosol  component.  Calcu- 
lations of  the  solar  aureole  for  540nm  and  for  large 
aerosol  optical  depth,  0.32  for  April  26,  1975,  show  an 
increase  in  aureole  intensity  of  0.1%  to  1.4%  for  scatter- 
ing angles  of  zero  and  fifteen  degrees,  while  for  a  small 
aerosol  optical  depth,  0.023  for  November  22,  1974,  an 
increase  of  0.6%  and  9.6%  for  the  same  scattering  angles 
was  observed.  The  correction  factor  is  less  important  at 
small  scattering  angles  because  the  strong  forward  peaking 
of  the  aerosol  scattering  makes  the  Rayleigh  component 
relatively  less  important.  As  can  be  seen  from  Table  5.3, 
the  correction  factor  is  largest  for  short  wavelengths, 
amounting  to  a  factor  of  two  at  312nm.  The  CDS  factor  is 
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meant  to  be  no  more  than  a  correction  factor,  a  small  cor- 
rection to  allow  the  single  scattering  calculation  to  more 
closely  model  reality,  so  at  ultraviolet  wavelengths  a 
true  multiple  scattering  calculation  should  probably  be 
made. 

5.4    Comparison  of  Measured  and  Calculated  Solar  Aureoles 

In  this  section  the  results  of  our  modified  single 
scattering  calculation  are  compared  in  a  number  of  graphs 
with  results  of  photographic  aureole  measurements.  Because 
the  measured  aureole  intensity  ratio  depends  on  an  accurate 
optical  density  measurement  of  the  neutral  density  four 
occulting  disc,  rather  a  high  density  for  precision,  we 
have  chosen  to  normalize  most  of  the  graphs  to  unity  at  one 
degree  scattering  angle  and  look  at  the  shape  of  the 
aureole  scattering  pattern. 

The  importance  of  the  index  of  refraction  is  shown  in 
Figure  5.8.  We   calculated  a  typical  aureole  pattern  for 
two  indices  of  refraction:  m  =  1.50+O.Oli,  which  is  a 
widely  used  value  for  the  average  index  of  refraction, 
and  1.33,  which  was  often  used  in  the  past.  Even  this  large 
a  change  in  index  of  refraction  has  little  effect  on  the 
shape  of  the  aureole,  although  on  an  absolute  scale  there 
is  a  substantial  change  in  magnitude,  partly  because  of 
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index  of  refraction  and  size  distribution  to 
the  solar  aureole. 


Figure  5.8 
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the  absorption  caused  by  the  O.Oli  imaginary  term.  In 
comparison,  the  effect  of  changing  the  size  distribution 
from  v=3.0  to  v=2.9  has  a  much  greater  effect. 

The  procedure  we  generally  followed  in  determining 
the  size  distribution  for  a  particular  day  was  first  to 
compare  experimentally  determined  aureoles  at  low  air  mass 
and  long  wavelength  (540nm  and  640nm)  to  the  corresponding 
single  scattering  calculations  for  an  array  of  v  and  a 
values.  The  amount  of  aerosol  used  for  the  calculation 
is  fixed  by  the  requirement  that  it  produce  the  same  aero- 
sol optical  depth  as  determined  by  the  solar  radiometer. 
Because  of  the  difficulty  in  obtaining  ozone  thicknesses 
from  the  Tallahassee  Dobson  station ,  we  used  the  average 
ozone  thickness  for  this  area  for  that  particular  month. 
In  Figure  5.9  we  show  how  strongly  variation  in  v  affects 
the  aureole.  There  was  no  difficulty  in  choosing  \»=3.2 
for  these  data.  Here  data  were  taken  on  April  26,  1975  in 
Gainesville,  which  is  about  sixty  miles  inland,  and  on 
Cedar  Key,  which  is  just  offshore  in  the  Gulf  of  Mexico, 
in  order  to  show  the  difference  between  an  inland  aerosol 
and  a  maritime  aerosol.  Unfortunately,  the  data  show  the 
aerosol  at  the  two  locations  to  be  identical,  even  though 
the  optical  depth  in  Gainesville  was  greater. 
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Comparison  of  measured  aureole  intensity  along 
a  horizontal  cut  through  the  sun  for  two  loca- 
tions with  single  scattering  calculations  for 
three  v  variations  in  the  size  distribution. 


Figure  5.9 
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The  calculations  for  Figure  5.9  were  made  with  a  size 
distribution  of  v=3.2  and  a=0.03  micron.  Because  the  sub- 
micron  particles  do  not  scatter  effectively,  little  differ- 
ence is  seen  whether  a=0.03  micron  or  a=0.10  micron  is 
used  in  the  calculation,  but  a  values  larger  than  this 
place  the  peak  of  the  size  distribution  into  the  effective 
light  scattering  domain  and  the  calculation  is  more  sensi- 
tive to  changes  in  a.  The  value  we  chose  is  in  agreement 
with  counter  measurements  which  indicate  that  the  size 
distribution  peaks  in  the  0.01  micron  to  0.10  micron  range 
usually.  Moreover,  if  the  size  distribution  peaks  at  a 
larger  radius,  the  allowable  ground  level  particle  density 
decreases  rapidly,  since  large  particles  are  very  effective 

attenuators;  yet  typical  particle  counter  measurements 

4 
in  Gainesville  show  about  10   particles  per  cubic  centi- 
meter and  an  average  background  mass  loading  of  about 
30  micrograms  per  cubic  meter.  The  size  distribution  used 
for  the  April  26  data,  which  was  a  day  of  large  aerosol 
optical  depth,  gives  a  ground  level  number  density  of 
4.4x10"*  particles  per  cubic  centimeter  and  75  micrograms 
per  cubic  meter,  niombers  consistent  with  Gainesville 
particle  counter  and  high  volume  sampler  measurements. 

In  Figure  5.10  we  plot  experimental  aureole  data  for 
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Figure  5.10 
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Comparison  of  measured  aureole  intensity  along 
a  horizontal  cut  through  the  sun  with  single 
scattering  calculation. 


Figure  5.11 
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Comparison  of  measured  aureole  intensity  along 
a  horizontal  cut  through  the  sun  with  single 
scattering  calculation. 


Figure  5.12 
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A  comparison  of  measured  aureole  intensity 
along  a  horizontal  cut  with  single  scattering 
calculation  for  two  albedo  values. 


Figure  5.13 


109 


November  22,  1974  at  four  wavelengths  on  an  absolute  scale, 
the  ordinate  being  the  ratio  of  aureole  to  direct  sxin 
intensity.  Note  the  increase  in  aureole  intensity  ratio 
with  decreasing  wavelength.  Of  course,  the  decreasing 
transmission  of  the  atmosphere  at  shorter  wavelengths  means 
that  the  actual  intensity  of  the  aureole  is  more  nearly 
constant  with  wavelength;  the  aureole  is  white,  not  blue. 

In  Figures  5.11,  5.12,  and  5.13  v/e  show  aureoles  for 
November  22  for  a  zenith  angle  of  54°  (air  mass  1.70)  at 
three  wavelengths.  Good  fits  to  the  data  for  this  date 
were  obtained  with  a  size  distribution  of  v=2.9  and 
a=0.03  micron.  The  importance  of  the  albedo  correction  is 
shown  in  Figure  5.13,  the  430nm  aureole.  The  calculations 
for  albedoes  of  0.0  and  0.25  do  not  differ  in  the  inner 
aureole,  but  at  ten  degrees  and  greater  the  fit  to  the 
data  is  significantly  better  for  an  albedo  of  0.25. 

Figures  5.14,  5.15,  and  5.16  show  the  solar  aureole 
for  the  same  date  and  wavelengths,  but  for  a  solar  zenith 
angle  of  75**  and  along  a  vertical  cut  below  the  sun  (toward 
the  horizon) .  The  two  sets  of  data  shown  for  the  540nm  and 
430nm  aureoles  are  in  each  case  from  different  negatives, 
one  made  with  a  50mm  lens  and  the  other  made  with  a  135rara 
lens.  The  calculated  aureole  predicts  an  increase  in  sky 
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Comparison  of  measured  aureole  intensity  along 
a  vertical  cut  below  the  sun  with  single  scat- 
tering calculation. 


Figure  5.14 
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Comparison  of  measured  aureole  intensity  along 
a  vertical  cut  below  the  sun  with  single  scat- 
tering calculation. 


Figure  5 . 15 
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Comparison  of  measured  aureole  intensity  along 
a  vertical  cut  below  the  sun  with  single  scat- 
tering calculation. 


Figure  5.16 
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brightness  near  the  horizon  at  640nm  and  540nin,  while  the 
data  shows  only  a  small  increase  at  the  horizon.  This 
discrepancy  could  be  caused  by  our  use  of  Eltermen's  model 
atmosphere.  These  aureoles  were  measured  on  an  exception- 
ally clear  day,  so  that  the  0-3  kilometer  aerosol  layer 
probably  did  not  contain  as  large  a  fraction  of  the  total 
aerosol  as  predicted  by  this  model,  and  the  low  level 
aerosols  are  responsible  for  the  near-horizon  brightness 
increase. 

Figures  5.17,  5.18,  and  5.19  are  sunset  aureoles, 
vertical  cuts  above  the  sun,  which  is  itself  only  two 
degrees  above  the  horizon.  The  photographs  here  were  not 
intended  for  vertical  cut  analysis,  so  for  the  640nm  and 
430nm  aureoles,  the  negatives  were  analyzed  along  a  vertical 
cut  adjacent  to  the  sun,  about  three  degrees  away  hori- 
zontally. Hopefully  this  cut  will  still  give  much  the  same 
aureole  pattern  except  for  the  first  few  degrees.  The 
agreement  between  data  and  calculation  was  rather  poor 
for  the  640nm  and  540nm  cases,  probably  because  of  the 
neglect  of  aerosol  multiple  scattering.  The  agreement  in 
the  430nm  aureole  is  better  because  Rayleigh  scattering  is 
relatively  more  important,  and  it  has  been  corrected  for 
multiple  scattering. 

Calculations;  for  the  309nm  aureole  are  shown  for  solar 
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to  sun  with  single  scattering  calculation. 


Figure  5.17 
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Figure  5.18 
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Comparison  of  measured  aureole  intensity  along 
a  vertical  cut  above  the  sun  and  adjacent  with 
single  scattering  calculation. 


Figure  5 . 19 
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zenith  angles  of  54*  and  75*  in  Figures  5.20  and  5.21.  As 
might  be  expected,  because  of  the  dominance  of  the  Rayleigh 
scattering  process  at  309nm,  the  single  scattering  calcu- 
lation including  the  CDS  correction  factor  does  a  much 
better  job  of  fitting  the  observed  solar  aureole  than 
does  the  single  scattering  calculation  shown  without  the 
correction  factor. 

For  comparison,  we  modified  the  Monte  Carlo  program 
written  by  Fuman  (8)  based  on  the  programs  of  Plass  and 
Kattawar  to  match  our  atmosphere  model  and  size  distri- 
bution. The  Monte  Carlo  program  follows  the  three  dimen- 
sional path  of  the  photon  as  it  penetrates  the  atmosphere, 
vmdergoing  scattering  and  absorption.  Channels  one  degree 
in  width  were  used  in  order  to  resolve  the  aureole;  unfor- 
tunately, use  of  such  a  large  number  of  channels  produces 
large  statistical  fluctuations  in  the  individual  channels, 
even  though  we  input  30,000  photons. 

The  Monte  Carlo  and  modified  single  scattering  calcu- 
lations agree  with  each  other  quite  well;  they  agree  with 
each  other  better  than  they  do  with  the  data.  Considering 
the  relative  costs,  the  modified  single  scattering  calcula- 
tion does  a  good  job  of  modeling  the  solar  aureole  even 
in  the  ultraviolet,  when  compared  with  a  Monte  Carlo 
calculation.  Only  for  the  sunset  calculations  at  longer 
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Comparison  of  measured  aureole  intensity  along 
a  horizontal  cut  with  single  scattering  calcu- 
lation, with  and  without  Rayleigh  multiple  scat- 
tering correction  factor,  and  a  Monte  Carlo 
multiple  scattering  calculation. 


Figure  5.20 
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Comparison  of  measured  aureole  intensity  along 
a  horizontal  cut  through  the  sun  with  single 
scattering  calculation. 


Figure  5.21 
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wavelengths  is  the  aerosol  multiple  scattering,  which 
we  neglect,  important  enough  to  cause  serious  error. 

5.5    Chromaticity  Diagrams 

The  color  of  the  sky  is  interesting  to  the  human 
observer,  and  by  using  a  system  of  coordinates  called  the 
tristimulus  values  a  quantitative  measure  of  the  color 
response  of  the  eye  to  the  calculated  spectral  sky  inten- 
sities may  be  derived.  The  standard  luminosity  functions 
x(X),  y(A)  ,  and  z'(X)  are  tabulated  by  Wyszecki  and  Stiles 
(46) ,  from  which  the  tristimulus  values  X,  Y,  and  Z  may 
be  calculated. 

X  =  /l(X)x(A)  dA  (5.43) 

Y  =  /l(X)y(A)  dA  (5.44) 

Z  =  /l(A)z(A)  dA  (5.45) 

From  these   the  chromaticity   coordinates  x,    y,    and  z   are 
obtained  through   the  expressions 

X  =  X/(X+Y+Z)  (5.46) 

y  =   Y/(X+Y+Z)  (5.47) 

z  =    Z/(X+y+Z)  (5.48) 
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Since  the  sum  of  the  three  coordinates  is  unity,  only  two 
coordinates  are  needed  to  specify  a  color. 

Figure  5.22  is  a  chromaticity  diagram  of  the  x  and 
y  coordinates.  The  center  of  the  diagram,  the  (0.333,0.333) 
point,  represents  white,  while  the  contours  outward  repre- 
sent increasing  saturation,  a  measure  of  the  degree  to 
which  a  pure  color  is  mixed  with  white  light,  reaching 
100%  saturation  on  the  outer  contour.  The  colors  are 
distributed  in  a  polar  direction  about  the  diagram  as 
indicated. 

Coordinates  for  two  different  zenith  angle  calcula- 
tions for  November  22  are  plotted  on  the  chromaticity 
diagram.  Direct  sunlight  for  the  54 *»  case  is  very  nearly 
white,  with  a  slight  yellow  tint  at  less  than  10%  satura- 
tion; while  the  88°  direct  sunlight  is  between  yellow  and 
orange  at  more  than  80%  saturation.  The  54°  aureole  is  of 
fairly  constant  hue,  blue,  but  going  from  20%  saturation 
at  0°  scattering  angle  to  more  than  30%  saturation  at  15°. 
The  sunset  aureole ,  corresponding  to  a  vertical  cut  through 
the  sky  up  from  the  horizon,  is  yellow  at  70%  saturation 
near  the  horizon,  decreasing  in  saturation  to  almost  white, 
then  becoming  slightly  blue  at  15°  above  the  sun.  It  is 
reassuring  that  the  calculated  sky  colors  at  least  behave 
in  a  reasonable  fashion. 


^ 


CHAPTER  6 


PHOTOGRAPHIC  ANALYSIS 


Photoelectric  detectors  are  most  widely  used  in 
studies  of  atmospheric  sky  brightness  because  the  photo- 
multiplier  tube  has  a  linear  response  over  a  large  dynamic 
range.  But  we  find  that  photographic  film  offers  several 
advantages  as  a  detector. 

The  major  advantage  is  the  rapidity  and  ease  with 
which  data  may  be  taken.  In  a  single  exposure  lasting  a 
thirtieth  of  a  second  we  record  sky  brightness  over  a  large 
area  of  sky,  an  area  26"  by  38"  for  a  50mm  lens,  at  very 
high  resolution.  The  modulation  transfer  function  for 
Panatomic-X  film  shows  that  it  will  resolve  a  spatial  fre- 
quency of  thirty  cycles  per  millimeter,  which  corresponds 
to  two  minutes  of  arc  resolution,  at  a  90%  response 
level.  And  if  one  accepts  a  30%  response,  the  resolution 
is  thirty  seconds  of  arc.  At  two  arc  minutes' resolution, 
one  frame  represents  a  900,000  point  array  stored  very 
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compactly  and  available  for  re-analysis  at  any  future  date. 
The  photographic  aureole  system  is  very  portable ,  indepen^ 
dent  of  any  power  supply,  so  that  aureole  data  may  be  taken 
from  a  variety  of  locations  with  a  minimum  of  preparation. 
Finally,  film  has  a  large  dynamic  range,  useful  over  three 
decades  of  intensity  and  linear  over  one  and  a  half  decades. 
Films  are  available  that  are  sensitive  from  about  110 Onm 
down  to  270nm  and  even  lower  for  low  gelatin  emulsions. 

6.1 Equipment 

We  photograph  the  skylight  near  the  sun  using  an  alt- 
azimuth mounted  35mm  camera.  The  mount  is  a  World  War  11 
searchlight  guidance  unit  which  may  be  finely  controlled 
in  azimuth  and  elevation.  The  multi-channel  solar  radio- 
meter is  usually  mounted  alongside  the  camera  to  simultan- 
eously measure  the  aerosol  optical  depth.  See  Figure  4.7. 

The  camera  we  use  is  a  Nikon  FTN,  which  is  excellent 
for  this  work  because  it  has  a  titanium  foil  shutter  rather 
than  the  usual  cloth  shutter  that  might  be  damaged  by  the 
direct  sunlight.  We  normally  use  a  50mm  lens  because  it 
covers  a  large  area  of  the  sky,  but  for  more  accurate  small 
angle  data  we  use  a  135mm  lens,  which  covers  only  10°  by 
15".  Either  lens  gives  accurate  information  to  within  one 
degree  of  the  sun's  center;  the  real  advantage  in  using  the 
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135inin  lens  is  that  it  spreads  the  sharp  decrease  in  aureole 
intensity  over  a  larger  area  of  film.  This  should  allow 
more  accurate  deteinnination  of  aureole  intensities,  but 
we  found  agreement  in  data  taken  with  the  two  lenses  to  be 
excellent  except  toward  the  edge  of  the  film,  where  insen- 
sitivity  of  the  scanning  rnicrodensitometer  to  slow  changes 
in  density  causes  discrepancy.  We  also  tested  a  SOOimn 
Nikon  reflex  (mirror)  lens  in  an  attempt  to  measure  the 
sub-one-degree  aureole,  but  the  long  hyperfocal  distance 
and  problems  with  internal  reflections  led  us  to  set  aside 
this  lens. 

The  ultraviolet  solar  aureole  at  309nm  we  photographed 
using  a  108mm  quartz  lens.  This  is  a  simple  one  element 
lens,  but  used  at  f/14  and  only  over  a  narrow  wavelength 
range,  its  performance  v/as  more  than  adequate. 

Direct  sunlight  is  blocked  from  entering  the  lens  by 
an  occulting  disc  held  about  one  meter  in  front  of  the 
camera.  This  is  necessary  because  the  sun  is  on  the  order 
of  ten  thousand  times  brighter  than  the  surrounding  aureole. 
The  larger  signal  would  easily  mask  the  smaller.  A  common 
solution  to  this  problem  is  use  of  a  modified  coronagraph, 
which  images  an  external  occulting  disc  onto  an  internal 
occulting  disc  to  remove  diffracted  light.  Rydgren  (47) 
has  measured  the  aureole  to  within  four  minutes  of  the  edge 
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of  the  sun  with  such  an  instrument.  In  using  only  an 
external  occulting  disc,  we  sacrifice  the  first  three- 
quarters  of  a  degree  of  the  aureole,  but  since  the  impor- 
tant features  of  the  aureole  occupy  the  one  to  ten  degree 
range,  there  is  no  compelling  reason  to  work  closer  to 
the  sun.  For  the  occulting  disc  we  use  a  neutral  density 
four  (ND4.0)  filter,  so  that  the  sun  appears  on  the  film 
attenuated  by  the  filter  along  with  the  aureole,  giving  us 
a  calibration  point  for  calculating  the  ratio  of  aureole 
intensity  to  direct  sunlight  intensity.  In  our  early  work 
we  used  a  Wratten  #96  ND  4.0  filter,  but  because  this 
material  increases  in  optical  density  at  shorter  wavelengths 
and  is  not  usable  at  309nm,  we  now  use  an  inconel  coated 
(reflection  type)  quartz  neutral  density  filter,  which  is 
spectrally  flat  from  250nm  to  2500nm  and  has  a  Callier  Q 
factor  of  1.005. 

Filters  are  used  to  limit  the  spectral  range  of  each 
aureole  photograph.  Interference  filters  centered  at  309nm, 
425nm,  and  850nm  were  each  less  than  lOnm  halfwidth,  and 
broadband  Wratten  absorption  filters,  #98,  #99,  and  #92 
centered  at  430nm,  540nm,  and  640nm  respectively  were  of 
about  BOnm  halfwidth.  The  #92  red  filter  is  actually  a 
long  pass  filter,  but  the  film  sensitivity  drops  to  zero 
at  690nm  to  give  an  effective  filter  cutoff  at  this  point. 
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The  filters  were  scanned  on  a  Gary  Spectrophotometer  to 
determine  filter  transmission  as  a  function  of  wavelength. 

We  normally  use  Panatomic-X  film  in  this  work  because 
of  its  low  speed,  fine  grain,  and  high  gamma.  The  film  is 
developed  in  D-76  at  68"?  (maintained  to  within  a  quarter 
of  a  degree)  for  eight  minutes.  This  film  is  sensitive  from 
below  300nm  to  690nm.  In  order  to  reduce  the  possibility 
of  long  wavelength  light  contributing  to  the  309nm  aureole, 
we  used  Spectrum  Analysis  No.  1  film  because  it  is  insen- 
sitive to  wavelengths  longer  than  470nm.  But  test  photo- 
graphs of  the  309nm  aureole  made  through  a  Mylar  filter, 
which  is  a  long  pass  filter  of  density  4.0  at  309nm  but 
only  density  0.5  at  320nm  and  above,  showed  that  leakage 
contributes  very  little  to  the  ultraviolet  aureole. 

Although  no  results  are  presented  here,  test  photo- 
graphs of  the  solar  aureole  at  850nm  were  made  using  High 
Speed  Infrared  film.  This  wavelength  is  in  a  water  vapor 
window,  ozone  no  longer  absorbs,  and  Rayleigh  scattering 
represents  an  optical  depth  of  only  0.02,  so  that  Mie 
scattering  remains  as  the  dominant  process .  In  some  ways 
tlHis  is  a  good  wavelength  at  which  to  study  aerosols,  but 
at  a  long  wavelength  the  small  particles  which  are  impor*r-' 
tant  for  visibility  in  the  atmosphere  are  no  longer 
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effective  scatterers,  and  the  aureole  is  no  longer  a  good 
indicator  for  particles  occupying  a  large  part  of  the 
size  distribution. 

6.2 Densitometric  Analysis 

It  is  essential  to  the  success  of  the  photographic 
aureole  study  to  be  able  to  accurately  relate  the  intensity 
of  a  point  in  the  aureole  to  the  density  of  a  corresponding 
point  on  the  film.  The  problem  is  in  two  parts:  the  geo- 
metric problem  of  imaging  the  sky  onto  a  flat  piece  of 
film,  and  the  sensitometric  problem,  which  will  be  dis- 
cussed in  the  next  section,  of  relating  the  exposure  the 
film  receives  to  the  subsequent  developed  density. 

A  simple  lens  (as  opposed  to  a  wide  angle  lens  for 
instance)  projects  light  entering  the  lens  from  the  direc-' 
tion  (6,(j))  to  a  point  on  the  film  in  the  direction  (tt+0, 
iT+(())  ;  i.e.  the  image  is  reversed  and  upside  down.  It  is 
natural  in  aureole  work  to  relate  positions  in  the  sky  to 
the  sun,  and  then  the  sun's  zenith  angle  and  azimuth  may 
be  calculated  for  a  given  time  and  location. 

Each  frame  is  scanned  on  a  Joyce  dual  beam  scanning 
microdensitometer  (Figure  6.1)  which  gives  density  and 
position  for  points  on  the  film.  An  x-y  coordinate  system 
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Joyce  dottble  beam  scanning  microdensitometer, 


Figure  6 . 1 
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is  established  with  zero  in  the  center  of  the  frame,  y 
positive  toward  the  zenith,  and  x  positive  toward  increasing 
azimuth.  The  coordinates  of  the  sun  are  x  ,y  and  a  point 
in  the  aureole  is  x  ,y  .  If  F  is  the  focal  length  of  the 

IT         hr 

lens,  the  zenith  angle  and  azimuth  angle  of  the  aureole 
point  are 


ZA  =  ZAg  +  arctan(yg/F)  -  arctan(y  /F)    (6.1) 


AZ   =  AZ  +  arctan(x  /F)  -  arctan(x„/F)    (6.2) 
p     s  p  s 


The  scattering  angle  is  given  by 


i|<  =  cos 


-1 


F-  +  XpX^  -f  y^y^ 

(F^+X^+y^^)l/2(p2^^2^y2)l/2 


(6.3) 


Light  received  at  the  film  plane  is  not  of  uniform 
irradiance,  but  varies  as  cos^G  where  Q   is  the  angle  rela- 
tive to  the  optical  axis.  At  points  off  the  axis  the  lens 
appears  foreshortened  into  an  ellipse,  so  that  its  apparant 
area  is  decreased  by  cos0.  Light  falling  on  the  film  at  an 
angle  is  spread  over  a  larger  area,  so  the  irradiance  on 
the  film  is  decreased  by  a  further  factor  of  cosG.  Finally, 
in  traveling  from  the  lens  to  the  film  the  irradiance 
decreases  as  1/r  and  since  r  is  inversely  proportional  to 
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cos9  this  accounts  for  the  final  factor  of  cos^9.  In  film 
coordinates,  the  angle  relative  to  the  optical  axis  is 


cos(e^)  =  J (6.4) 

S        /  ,,2  ,..2  ...2 


s  ■'s 


cos(e  )  — — (6.5) 

•^  F^+xf^+y!, 
P  P 


6.  3 Sensitometry 

Sensitometry  relates  the  exposure  received  by  the  film 
to  its  developed  density.  The  equation  relating  exposure 
to  blackening  of  film,  known  as  the  characteristic  curve, 
was  first  carefully  studied  in  1890  by  Hurter  and  Driffield 
(48) ,  and  the  characteristic  curve  is  frequently  referred 
to  as  the  H  and  D  curve. 

In  order  to  determine  the  characteristic  curve  for  a 
particular  film,  developer,  and  light  source  combination, 
it  is  necessary  to  deliver  to  the  film  a  series  of  graded 
exposures  covering  the  range  of  sensitivity  of  the  film. 
Since  exposure  is  the  product  of  irradiance  multiplied  by 
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the  duration  of  exposure,  the  simplest  method  of  doing  this 
would  be  to  expose  at  constant  irradiance  for  varying  times. 
But  because  of  reciprocity  failure,  a  high  intensity  short 
duration  exposure  does  not  produce  the  same  density  as  a  low 
intensity  long  duration  exposure.  A  good  technique  for  pro- 
ducing time  variation  characteristic  curves  is  to  use  a 
rotating  stepped  sector  wheel.  At  low  rotation  rates  the 
intermittency  effect  makes  this  unreliable,  but  above  a 
certain  rotation  rate  called  the  fusion  frequency  this  is 
a  useful  method. 

The  alternative  to  time  variation  is  intensity  (irrad- 
iance) variation.  A  simple  form  of  intensity  variation  is 
exposure  of  film  through  a  calibrated  step  wedge.  Because 
we  require  a  sensitometer  that  is  useful  over  a  wide  range 
of  wavelengths,  v/e  use  a  tube  sensitometer  similar  to  the 
Harvard-Smithsonian  tube  sensitometer  described  by  Latham 
(49).  In  this  device  a  plate  is  suspended  in  an  integrating 
sphere  to  give  a  very  evenly  illuminated  surface  which  is 
viewed  by  a  tube  stack,  each  tube  acting  as  a  pin-hole 
camera  giving  an  irradiance  at  the  front  plate  proportional 
to  the  area  of  the  entrance  aperture  of  that  txibe.  (See 
Figure  6.2.)  Latham  shows  that  the  stack  geometry  deter- 
mines the  variation  in  intensity  across  the  exit  aperture. 
The  fractional  variation  from  center  to  edge  is  approxi- 
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Tube  sensitoineter,  "mini-sensitometer" 
version,  showing  integrating  chamber, 
reflection  plate,  and  tube  stack  with 
varying  aperture  front  plate. 


Figure  6.2 
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mately  given  by 


^-     =  2(R/t)2  (6.6) 


where  R  is  the  radius  of  the  spot  and  t  is  the  height  of 
the  tube  stack.  This  gives  a  fractional  variation  for  the 
largest  tube  in  our  sensitometer  of  0.006,  about  half  a 
percent. 

The  integrating  "sphere"  in  our  sensitometer,  shown 
in  Figure  6.3,  is  rectangular,  30cinx45cmx50cin.  The  tube 
stack  is  designed  to  give  a  strip  of  35iTun  film  twenty-two 
exposures  in  two  rows  of  eleven,  covering  an  exposure 
range  of  three  orders  of  magnitude  in  approximately  equal 
logarithmic  steps.  Each  spot  is  a  quarter  inch  in  diameter 
and  well  separated  to  avoid  overlap  due  to  halation.  The 
integrating  chamber  is  coated  with  magnesium  oxide,  a  good 
diffuse  reflector  over  a  wide  range  of  wavelengths. 

The  sensitometer  does  not  have  an  integral  light 
source.  The  entrance  aperture  has  a  manual  slide  for  use  as 
a  shutter  and  can  be  equipped  with  a  leaf  shutter  for  very 
short  exposures.  At  various  times  we  have  used  a  500  watt 
lamp,  a  1000  watt  quartz-iodine  lamp,  a  high  intensity 
Bausch  and  Lomb  monochromator  for  studying  spectral  varia- 
tion in  gamma  in  the  ultraviolet,  and  natural  sunlight 
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Tube  sensitometers  used  to  obtain  H  and  D  curves 
for  35inm  film. 


Figure  6.3 
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concentrated  by  a  ten  inch  parabolic  mirror  for  very  short 
exposure  times . 

In  addition  to  our  main  sensitometer  we  built  a  "mini- 
sensitometer" ,  also  shown  in  Figure  6.3,  only  ten  centi- 
meters in  diameter  and  with  a  nine  hole  tube  stack  to  study 
the  response  of  human  skin  to  ultraviolet.  The  throughput 
of  the  small  sensitometer  is  almost  an  order  of  magnitude 
greater  than  the  large  version. 

The  sensitometer  was  calibrated  using  a  photomulti- 
plier  detector,  and  it  was  found  that  the  agreement  between 
the  aperture  area  and  the  photomultiplier  readings  was 
generally  good.  As  an  additional  check  film  was  exposed 
so  that  the  linear  part  of  the  HD  curve  fell  on  the  differ- 
ent parts  of  the  exposure  range;  any  spots  not  properly 
calibrated  would  consistently  deviate  from  the  straight 
line. 

6.4    Analytic  II  and  D  curves 

The  first  attempt  to  quantitatively  relate  the  trans- 
parency of  an  image  to  its  exposure  was  made  by  Abney  (50) 
who  suggested  that  the  relation  was 

T  =  3-u(log{E/i})^  (6^7) 
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where  the  transmittance  T  is  related  to  the  density  D 
and  the  opacitance  0  by 

D  =  log(l/T)  =  log(O)  (6.8) 

Hurter  and  Driffield  developed  the  experimental  methods  for 
determining  the  characteristics  of  a  photographic  layer  and 
the  nomenclature  for  its  optical  properties  which  has  since 
been  generally  adopted.  They  foiind  that  at  low  exposure 
levels  the  sensitive  layer  builds  up  density  nearly  linear- 
ly with  exposure  until  the  main  part  of  the  curve  is 
reached,  where  density  increases  as  the  logarithm  of  the 
exposure.  At  high  exposure  levels  nearly  all  the  silver  is 
deposited  and  a  saturation  is  reached  where  density 
increases  more  and  more  slowly  with  exposure.  This  is  the 
familiar  sigmoid  curve  as  in  Figure  6.4  that  is  obtained 
when  density  is  plotted  against  the  logarithm  of  the  expo- 
sure. Over  the  main  part  of  the  curve  the  relation  is 
linear  of  slope  y ,    and  the  exposure  corresponding  to  the 
extrapolation  of  the  straight  line  to  the  base  density  is 
the  inertia  of  the  film,  labeled  a  in  Figure  6.4.  For  the 
straight  line  part  of  the  HD  curve  they  found  that 

D  =  YdogE  -  logi)  (6.9) 
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where  i,  the  inertia,  is  expressed  in  exposure  units. 

The  equation  of  Hurter  and  Driffield  holds  only  for 
the  central  part  of  the  curve  and  has  no  obvious  fundamen- 
tal physical  basis.  A  number  of  workers  since  have  attempted 
to  deduce  an  expression  for  the  entire  curve,  among  them 
Slade  and  Higson  (51)  and  Silberstein  and  Trivelli  (52). 
Such  basic  approaches  met  with  varying  degrees  of  success, 
but  at  the  price  of  complexity  in  the  analytic  expressions. 
So  when  Vaucouleurs  (53)  recently  pointed  out  a  simple 
analytic  expression,  originally  discovered  by  Baker,  for 
the  linearization  of  the  toe  of  the  characteristic  curve, 
his  paper  met  with  a  great  deal  of  interest. 

In  terms  of  the  opacitance  to 


E  =  EqJ^  (6.10) 


where 


0)  =  ^  -  1  =  10°"°L  -   1  (6.11) 


D-  corresponds  to  the  base  density  of  the  film  and  N  to 
the  inverse  of  the  slope,  l/y.  A  variant  on  this  form  is 

given  by  Mees  (54) 

D-Dl 
E  =  Eq(1Q  u   -  i)N  (6.12) 
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where  u  is  given  by  the  density  of  the  inertia  point  by 
u  =  D./2.  Green  and  McPeters  (55)  have  studied  a  further 
generalization  that  linearizes  the  shoulder  of  the  charac- 
teristic curve  as  well  as  the  toe,  by  defining  a  quantity 
similar  to  the  opacitance  for  the  saturation  density 


(0=1-  10^<°-°u>  (6.13) 

u 


where  D  is  the  saturation  density  and  3  is  a  measure  of 
the  asynunetry  of  the  curve.  If  we  define  a  total  opacitance 

fi  =  (o/Uy  (6.14) 

then  the  entire  characteristic  curve  is  given  by 


E  =  Eq  fi^  (6.15) 


or,  in  terms  of  density 


E  =  Eq{(10°-'^L  -  1)/(1  -  10^(°-°u>)}N        ^g  ;Lgj 


Characteristic  curves  for  Spectrum  Analysis  No.  1  and 
Panatomic-X  film  are  shown  in  Figure  6.4  along  with  fits 
given  by  this  equation. 
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Characteristic  curves  for  Panatomic-X  at  four  wave- 
lengths are  shown  in  Figure  6.5.  These  curves  were  used 
with  the  April  26,  1975  aureole  data.  The  film  was  exposed 
in  the  mini-sens itometer  to  sunlight  through  the  same 
filters  used  in  shooting  the  aureoles  and  for  the  same 
exposure  times.  The  film  from  the  sensitometer  was  then 
processed  in  the  same  tank  with  the  aureole  film.  We  found 
that  we  could  obtain  good  fits  to  this  data  with  6  set  to 
unity.  The  best  fit  parameters  for  the  four  wavelengths 
are  given  in  Table  6.1. 

The  most  general  form  of  the  equation  is 

^  a(D-DL)^  _  i}N 

E  =  E  -i^ fl-  (6.17) 

{1  -  ioe(D-Du)l^jM 

where  the  parameters  a,3/N,M,v,  and  jj  may  be  set  to  other 
than  unity  as  required. 

As  might  be  indicated  by  the  procedure  described  for 
the  April  26  data,  for  greatest  accuracy  it  is  necessary  to 
match  sensitometric  procedures  to  the  film  use  as  closely 
as  possible.  Reciprocity  failure  requires  exposure  time 
matching,  and  because  an  HD  curve  can  change  sharply  with 
wavelength,  the  spectral  quality  of  the  light  source  is 
important.  Naturally  the  temperature  and  freshness  of  the 
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TABLE  6.1 


Best  fit  parameters  for  HD  curves  of  Figure  6.5. 


parameters 

309 

430 

540 

640 

DL 

0.245 

0.245 

0.245 

0.245 

DU 

3.00 

3.04 

2.96 

2.96 

N 

0.759 

0.649 

0.653 

0.613 

=0 

52.9 

8.83 

5.75 

5.49 
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developer,  the  time  of  development,  and  the  agitation 
procedure  are  all  important.  Finally,  the  same  densito- 
metric  procedure  must  be  used  on  all  negatives  and  sensi- 
tometer  strips.  Because  the  Callier  Q  factor  for  film 
varies  greatly  with  density  and  grain  size  and  may  be  1.4 
or  higher  for  middle  densities,  two  densitometers  that 
admit  differing  amounts  of  scattered  light  will  indicate 
different  densities. 

Despite  all  precautions ,  we  found  a  certain  amount 
of  scatter  in  our  characteristic  curve  data.  Especially 
in  the  high  density  spots,  we  sometimes  noted  density  dif- 
ferences from  one  edge  of  the  spot  to  the  other  which  we 
attribute  to  non-uniformities  in  agitation.  We  develop 
the  film  in  a  two  reel  stainless  steel  tank,  using  the 
standard  procedure  of  agitation  at  thirty  second  intervals, 
To  achieve  better  uniformity,  a  better  system  of  agitation 
such  as  nitrogen  burst  agitation  or  a  mechanical  agitator 
should  be  used.  But  we  feel  that  the  errors  involved  are 
not  large  since  a  smooth  curve  will  pass  through  enough 
points  to  compensate  for  individual  errors. 
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6.  4 Conclusions 

We  have  demonstrated  in  this  dissertation  the  useful- 
ness of  a  nuinber  of  techniques  for  computations  of  atmos- 
pheric light  scattering. 

The  use  of  analytic  rather  than  numerical  models  is 
a  strong  technique  that  should  be  more  widely  used  in 
atmospheric  optics  problems.  As  we  have  shown  a  few  rather 
simple  analytic  forms  can  be  applied  quite  flexibly  to 
reasonably  smooth  distributions;  only  highly  structured 
distributions,  such  as  the  extraterrestrial  solar  flux  or 
the  spectral  ozone  absorption  coefficients,  require 
numerical  modeling. 

We  were  able  to  calculate  aerosol  optical  depths  from 
filter  radiometer  readings  in  the  middle  ultraviolet  by 
fitting  an  analytic  function  to  the  calculated  "effective" 
optical  depth  seen  by  the  instrument.  Our  technique  makes 
possible  routine  measurement  and  calculation  of  the  aerosol 
optical  depth  in  the  ultraviolet,  something  previously  not 
normally  done  with  filter  radiometers. 

We  explored  the  use  of  analytic  functional  forms  for 
the  Mie  scattering  phase  functions  for  various  size  distri- 
butions, but  further  work  is  necessary  for  this  to  be 
generally  useful.  Some  form  of  eigenvector  technique  to 


146 


generate  the  phase  function  for  an  arbitrary  size  distri- 
bution is  needed. 

Many  others  have  studied  the  solar  aureole,  but  our 
technique  of  measuring  the  sky  to  sun  intensity  ratio 
through  use  of  a  neutral  density  occulting  disc  is,  we 
believe,  unique.  The  sunset  sky  brightness  has  recently 
been  calculated  for  realistic  atmosphere  models.  We  concen- 
trate on  the  near-sun  sky  brightness  and  attempt  to  match 
actual  sunset  aureoles  to  calculations  that  use  an  atmos- 
pheric model  that  accurately  predicts  the  observed  aureoles 
at  smaller  zenith  angles  on  a  given  day.  That  we  correctly 
calculate  some  features  of  the  sunset  aureole  is  promising, 
but  aerosol  multiple  scattering  apparcintly  must  be  in- 
cluded to  correctly  match  the  actual  sunset  aureole. 

The  modified  single  scattering  calculation  which  we 
have  developed  is  sufficiently  flexible  to  match  a  wide 
range  of  atmospheric  conditions  and  in  many  cases  offers 
significant  improvement  over  a  simple  single  scattering 
calculation. 
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